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ABSTRACT 


A  study  was  conducted  to  evaluate  the  accuracy  of  a  commercial  computational  fluid 
dynamics  (CFD)  code  (CFDRC-ACE+)  for  predicting  incompressible  air  jet  flows  with 
simple  geometries.  Specifically,  the  axis-symmetric  and  two-dimensional  heated  air-jets  were 
simulated  using  a  standard  k-z  turbulence  model.  These  CFD  predictions  were  directly 
compared  to  an  extensive  compilation  of  experimental  data  from  archive  literature.  The  round 
jet  results  indicated  that  the  code  over-predicted  the  velocity-spreading  rate  by  24%  and  the 
temperature-spreading  rate  by  29%.  In  addition,  the  centerline  velocity  and  temperature  decay 
rates  were  also  over-predicted  by  21%  and  30%,  respectively.  The  geometric  and  kinematic 
virtual  origins  were  over-predicted,  as  well,  by  approximately  7.5  diameters  for  the  velocity 
profiles  and  10.5  diameters  for  the  temperature  profiles.  The  planar  jet  simulation  was 
generally  closer  to  experimental  data  ranges,  with  an  under-prediction  of  the  velocity¬ 
spreading  rate  of  approximately  17%  with  an  over-predicted  temperature-spreading  rate  of 
12%.  The  centerline  velocity  and  temperature  decay  rates  were  both  under-predicted  at  22% 
and  27%,  respectively.  Again,  the  geometric  and  kinematic  virtual  origins  were  over¬ 
predicted  by  approximately  7.5  slot  heights  for  the  velocity  profiles  and  10.5  slot  heights  for 
the  temperature  profiles. 
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INTRODUCTION 


A.  MOTIVATION 

The  characteristics  of  jet  development  are  very  important  in  many  military  and  industrial 
applications  including  eductors,  thrast  augmenting  ejectors  and  exhaust  suppression  units.  For 
example,  current  efforts  are  underway  to  improve  present  shipboard  exhaust  eductors  without 
adversely  effecting  other  parts  of  the  ship.  Former  redesign  of  these  types  of  systems  has  mainly 
been  accomplished  through  extensive  scale  model  testing  of  each  new  design  proposal.  The 
desire  to  reduce  developmental  costs  of  design  test  configurations  has  stimulated  an  increased 
interest  in  Computational  Fluid  Dynamics  (CFD)  design  studies.  By  using  CFD  analysis,  only 
the  most  favorable  design  modifications  warrant  further  experimental  testing. 

B.  CHALLENGE 

CFD  simulations  help  researchers  to  understand  and  visualize  the  complex  flow  patterns 
before  costly  experiments  are  performed.  However,  CFD  codes  are  typically  challenged  when 
predicting  flows  that  are  dominated  by  turbulent  shear  mixing.  The  general  conservation 
equations  of  continuity,  momentum  and  energy  can  be  solved  directly  by  using  Direct  Numerical 
Simulation  (DNS)  on  supercomputers.  However,  the  extremely  small  grid-cell  size  and  time 
steps  required  to  resolve  the  smallest  scales  of  turbulence  (Kolmogorov  scale)  limit  DNS 
calculations  to  low/moderate  Reynolds  number  (Re)  flows  only. 
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In  order  to  reduce  the  amount  of  computational  resources  needed,  turbulence  models 
have  been  developed  to  mimic  the  effects  of  turbulence  (Reynolds  stresses)  that  are  contained 
within  the  conservation  equations.  This  approach  permits  the  use  of  smaller  computers  to  model 
complex  and  high-speed  flow  configurations.  However,  turbulence  models  merely  approximate 
the  Reynolds  stresses  and  are  usually  limited  to  a  small  number  of  simple  flow  regimes.  Often 
‘turbulence  model  constants’  must  be  adjusted  to  make  the  model  accurately  predict  other  types 
of  flows. 

Depending  on  the  turbulence  model  used,  any  number  of  empirically  derived  constants 
may  be  required  to  properly  model  the  flow.  Overall,  the  process  of  correctly  modeling  the 
Reynolds  stresses  of  a  flow  configuration  is  extremely  difficult  and  problematic.  Ideally,  a 
calculation  derived  from  the  flow  geometry  itself  would  be  the  most  desirable  option  when 
calibrating  a  turbulence  model.  However,  only  empirically  derived  constants  from  relatively 
simple  experiments  have  been  able  to  provide  such  inputs.  Unfortunately,  the  highly  complex 
flow  patterns  that  typically  arise  in  nature  have  precluded  these  experimental  constants  from 
being  truly  universal.  Therefore,  a  “trial  and  error”  methodology  is  often  employed  to  adjust  the 
turbulence  model  constants  for  reasonably  accurate  predictions  in  complex  flows. 

Commercial  codes,  designed  to  model  many  types  of  flows  with  average  operator  skill, 
often  make  certain  concessions  to  provide  acceptable  overall  predictions.  Frequently,  the 
constants  are  an  integral  part  of  a  commercial  code  and  can  not  be  changed  to  fit  the  flow  regime 
under  consideration.  On  the  other  hand,  ‘wall  functions’  and  other  additional  turbulence  features 
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are  commonly  used  to  modify  or  correct  the  local  flow  conditions  as  necessary.  Overall,  the 
process  of  using  turbulence  models  is  highly  challenging  and  the  results  are  sometimes  less  than 
desired. 

C.  OBJECTIVE 

To  eventually  predict  complex  multiple-jet  type  flows,  a  single  free  jet  must  first  be 
correctly  modeled  without  the  interference  of  walls  or  adjacent  jets.  Accurate  simulation  of  the 
simple  ‘free  jet’  is  crucial  for  the  continued  use  of  the  CFD  code  and  turbulence  model.  With 
this  baseline  of  code  experience  and  knowledge,  further  research  can  then  progress  towards 
developing  and  configuring  more  complex  flow  configurations. 

The  objective  of  this  research  is  to  determine  the  performance  and  sensitivity  of  CFD 
Research  Corporation’s  numerical  analysis  code  (CFD-ACE+)  to  predict  simple  nearly 
incompressible  free  turbulent  jet  flows.  In  particular,  heated  air  axisymmetric  (round)  and  two- 
dimensional  (planar)  jets  were  modeled  and  compared  to  tabulated  experimental  data.  In 
addition,  a  critical  review  of  experimental  data  available  in  the  archive  literature  is  also  provided. 
The  difficulty  in  predicting  these  simple  flows  has  proven  to  be  one  of  the  more  challenging 
problems  faced  by  many  turbulence  modelers  over  the  last  five  decades.  The  CFD-ACE+  code 
experience  and  jet  behavior  knowledge  provided  by  this  research  was  meant  to  create  a  solid 
foundation  for  further  investigations. 
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n.  BACKGROUND  ON  FREE  JETS 


A.  PREFACE 

The  turbulent  jet  has  been  the  subject  of  many  extensive  studies,  both  theoretical  (i.e., 
Abramovich  1963,  Hinze  1975,  Townsend  1976,  Schlichting  1979,  Schetz  1980,  List  /  Rodi 
1982),  and  experimental  (i.e.,  Hinze  /  van  der  Hegge  Zijnen  1949,  Wygnanski  /  Fiedler  1969, 
Rodi  1975,  Panchapakesan  /  Lumley  1993,  Hussein  et  al  1994).  Being  both  basic  and  important 
in  many  applications,  ‘simple  geometry’  jets  are  actually  quite  complex  and  multifaceted.  From 
the  earliest  investigations,  the  turbulent  jet  has  been  dissected  into  different  discrete  regions  for 
in-depth  analysis.  A  brief  description  will  be  given  here  to  familiarize  the  reader  with  this  type 
of  flow. 

A  free  jet  is  formed  when  a  fluid  is  discharged  from  an  orifice  or  nozzle  into  quiescent, 
co-flowing,  or  counter-flowing  environments.  Co-flows  and  counter-flows  are  not  considered  in 
this  investigation.  By  analyzing  only  subsonic  flows  with  Mach  Numbers  less  than  0.3, 
compressibility  effects  can  be  ignored.  In  addition,  only  jets  with  top-hat  velocity  / 
temperature  exit  profiles,  negligible  buoyancy  effects,  and  stagnant  surroundings  are  considered. 
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As  shown  in  Figure  1,  the  jet  is  usually  divided  into  three  main  regions.  Each  region  has 
its  own  unique  characteristics  that  set  it  apart  from  the  others.  The  virtual  origin  is 
conceptualized  as  an  “ideal”  jet  source  that  produces  the  same  self-similar  jet  effects  without  the 
transient  regions  (potential  core  and  intermediate).  The  “geometric”  virtual  origin,  as  indicated 
in  Figure  1,  is  obtained  by  extending  the  linear  spreading  profile  of  the  jet  backward  to  the  jet’s 
centerline.  The  distance  from  the  jet  center  (where  the  velocity  is  Uc)  to  the  half-velocity  point 
(Uc/2)  usually  defines  the  jet’s  spreading  profile  length.  The  “kinematic”  virtual  origin,  on  the 
other  hand,  is  acquired  from  the  inverted  centerline  velocity  decay  plot  (Ue/Uc).  Again,  a  line  is 
extended  backward  from  the  linear  self-similar  decay  region  to  the  X-axis  to  indicate  the  virtual 
origin  as  shown  in  Figure  2. 
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Centerline  Velocity  Decay 


Figure  2.  Inverted  Centerline  Velocity  Decay  Plot 
B.  POTENTIAL  CORE  REGION 

The  potential  core  region  (Figure  1),  closest  to  the  orifice,  is  highly  influenced  by  the 
outlet  conditions  of  the  orifice  as  well  as  the  surrounding  medium.  The  centerline  streamwise 
velocity  remains  approximately  constant  in  this  region;  it  may  however  exhibit  a  slight  vena 
contracta  at  higher  flow  rates.  The  significant  shear  distortions  in  the  flow  field  rapidly  generate 
instabilities  and  create  turbulent  energy  in  the  shear  layers  and  greatly  effect  the  initial  spreading 
rate  of  the  jet.  Thought  by  many  early  investigators  to  be  totally  incoherent  fluid  motions,  this 
turbulent  flow  has  been  proven  to  have  ‘organized-structures’  (coherent  structure)  superimposed 
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on  its  chaotic  (or  random)  motion.  This  large-scale  turbulence  generated  at  the  jet  boundaries  is 
most  responsible  for  noise  and  initial  mass  entrainment.  List/Rodi  (1982)  review  several  modes 
of  wave/vortex  growth  and  instabilities  and  their  descriptions  will  not  be  repeated  here. 

Generally,  by  lowering  inlet  core  turbulent  intensity  (Tu)  and  flattening  velocity  profiles 
(top  hat),  longer  potential  core  lengths  are  realized  in  this  region.  Also  increasing  the  exit 
velocity  (effectively  lowering  Tu  and  flattening  the  velocity  profile)  can  enlarge  the  potential 
core  length. 

C.  INTERMEDIATE  REGION 

As  entrainment  develops  on  its  periphery,  the  potential  core  disappears  into  the 
intermediate  (transitional)  region.  The  width  of  this  transitional  region  spreads  rapidly  outward 
as  the  jet  progresses  down  stream  and  the  average  centerline  velocity  decreases.  This  re-shaping 
of  the  mean  flow  tends  to  force  the  jet  into  self-similarity  (self-preservation). 

D.  SELF-SIMILAR  REGION 

Eventually,  the  flow  attains  self-similarity  in  which  the  various  terms  of  the  momentum 
equation  (convection,  diffusion,  production,  etc.)  maintain  constant  ratios  so  that  the  various 
processes  are  in  dynamic  equilibrium,  each  one  changing  downstream  at  the  same  rate  as  the 
others.  In  this  region,  dimensionless  velocity  and  length  scales  (a  function  of  one  geometrical 
variable  only)  can  completely  describe  the  flow.  The  flow  is  also  said  to  be  independent  of  the 
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initial  jet  conditions  including  orifice  shape  and  turbulence  intensity.  The  width  of  this  region 
seems  to  grow  at  a  linear  rate  as  the  flow  diminishes.  Mean  velocity  profile  similarity  is  usually 
reached  within  about  20-30  jet  diameters  downstream  while  turbulence  velocity  fluctuations 
don’t  obtain  “true  similarity”  until  70  diameters  or  greater  depending  on  the  type  of  flow.  Many 
books,  such  as  Abramovich  (1963),  Launder  (1975)  and  others,  discuss  similarity  arguments  in 
much  greater  detail. 

Jets  can  obtain  “exact”  similarity  only  if  they  issue  into  stagnant  or  quiescent 
surroundings  (as  opposed  to  co-flow  or  counter-flow)  and,  in  general,  are  theorized  to  have  the 
following  approximate  streamwise  characteristics  under  the  assumptions  of:  steady  flow, 
negligible  body  force  and  negligible  pressure  gradient  (uniform  static  pressure). 

For  axisymmetric  round  jets,  the  magnitude  of  axial  velocity  at  the  center  of  the  jet  (Uc) 
varies  as  the  inverse  of  the  axial  distance  (X)  from  the  jet  exit  nozzle.  On  the  other  hand,  the 
radius  at  which  the  axial  velocity  of  the  jet  profile  is  half  of  the  centerline  value  (R1/2)  varies 
directly  with  X. 

Centerline  Velocity  Decay:  Uc  =  X‘*  (1) 

Half-Velocity  Radius:  R1/2 « X  (2) 
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Planar  jets  exhibit  slightly  different  characteristics  due  to  their  two-dimensional  nature. 
The  centerline  velocity  (Uc)  and  transverse  half- velocity  length  (Y1/2)  vary  in  the  axial  direction 
(X)  as  follows. 


Centerline  Velocity  Decay: 

Uc  « 

(3) 

Half- Velocity  Length: 

Y,/2-X 

(4) 

E.  TEMPERATURE  PROFILE 

The  turbulent  mixing  within  a  jet  causes  a  transfer  of  fluid  particles  in  a  direction  normal 
to  the  main  flow.  The  process  is  complex  and  involves  momentum  flowing  away  from  the  main 
stream  while  the  surrounding  fluid  is  being  entrained  toward  it.  In  addition,  there  can  be  a 
transfer  of  heat  or  diffusion  of  the  temperature  field  into  the  surroundings.  As  shown  in  Figure 
3,  the  temperature  profile  is  considerably  broader  than  the  velocity  profile  due  to  the  higher 
transverse  transport  of  thermal  energy  by  turbulence  than  momentum  for  gas  jets.  The  scalar 
values  of  temperature,  internal  energy  and  species  concentration  all  behave  in  the  same  way  for 
these  types  of  flow. 

Abramovich  (1963),  Hinze  (1975),  Schlichting  (1979)  and  others  describe  heat  transfer 
theories  within  the  jet  that  are  based  on  several  simplifying  assumptions.  Although  in  general, 
the  exact  jet  characteristics  are  not  presented  in  the  texts  and  the  authors  apply  their  (or  other’s) 
experimental  data  to  obtain  the  required  empirical  constants.  However,  each  jet  tends  to  be 
unique  and  the  limited  data  used  by  each  author  are  only  applicable  for  that  specific  jet  which 
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may  only  be  similar  to  one  or  two  other  experiments.  Even  though  these  limited  data  (usually 
dated)  will  support  the  author’s  general  conclusions,  many  other  experiments  with  comparable 
setups  do  not  obtain  the  same  results.  Obviously  not  all  factors  (like  coherent  structures,  etc) 
have  been  properly  controlled  in  the  experiments  and  incorporated  into  equations  (1  to  4)  and 
further  investigation  is  certainly  warranted.  As  will  be  discussed  in  Chapter  HI,  factors  like 
density  ratio  and  facility-related  disturbances  can  greatly  effect  the  developing  characteristics  of 
free  jets. 


■n=R(X-Xo) 


Figure  3.  Typical  Velocity  and  Temperature  Profiles 
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in.  CRITICAL  REVIEW  OF  EXPERIMENTAL  RESULTS 


A.  OVERVIEW 

It  is  essential  to  have  reliable  experimental  data  against  which  CFD  predictions  can  be 
compared.  Free  turbulent  jets  have  been  measured  fairly  extensively  throughout  the  years  and 
plenty  of  data  are  available.  The  purpose  of  this  chapter  was  to  review  these  data,  check  its 
reliability,  and  select  target  values  for  a  comparison  with  the  predictions  from  the  CFD  code.  As 
the  reader  will  soon  discover,  this  task  is  not  as  simple  as  it  first  may  appear. 

B.  MEASUREMENT  TECHNIQUES 

In  order  to  make  judgements  on  the  particular  quality  of  the  measurements;  we  will 
briefly  examine  the  merits  and  shortcomings  of  various  measurement  techniques  used  to  obtain 
the  referenced  data.  In  most  free  jet  flows,  the  static  pressure  varies  very  little  across  the  flow. 
Therefore,  measurement  of  total  pressure  by  a  Pitot  tube  (or  similar)  can  be  sufficient  to 
determine  the  mean  velocity  (U).  The  Pitot  tube,  due  to  its  slow  response  however,  can’t 
measure  the  fluctuating  velocities  that  represent  the  most  serious  source  of  error.  When 
turbulence  fluctuations  are  high  (>20%  -  like  in  the  shear  layer),  the  instantaneous  yaw  angles 
can  be  large  enough  to  affect  the  probe  reading  without  any  reliable  correction  available. 
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The  hot-wire  anemometer  has  a  sufficiently  fast  response  to  follow  the  turbulence 
fluctuations  and  has  become  the  one  of  the  more  popular  instruments  for  measuring  turbulent 
flows.  However,  many  sources  of  error  are  associated  with  the  hot-wire  measurement  technique 
and  not  all  early  researchers  have  properly  accounted  for  them  as  described  by  Hussein  et  al 
(1994)  and  Panchapakesan  /  Lumley  (1993).  The  more  common  errors  include  improper  unit 
calibration,  turbulence  interference  between  adjacent  wires  (multiple  wire  configurations),  and 
failure  to  realize  that  a  significant  fraction  of  reverse  flow  exists  in  the  outer  shear  layers  of  jets. 
These  problems  have  introduced  many  uncertainties  and  the  general  reliability  of  the  data 
essentially  decreases  with  increasing  turbulence  level. 

The  relatively  new  laser  technologies  offer  much  promise  for  high  accuracy  data.  By 
applying  a  known  fi-equency  shift  to  the  laser  beams,  an  offset  frequency  is  created  that  is 
increased  or  decreased  by  the  Doppler-shifl  phenomena.  Laser— Doppler  Anemometers  (LDA) 
are  one  of  the  very  few  velocimeters  capable  of  measuring  flow  reversals.  However,  concerns 
still  arise  about  proper  seeding  material  size  and  dispersion  required  to  truly  represent  the  flow. 
Presently,  only  very  few  jet  flow  configurations  have  been  re-examined  using  this  state-of  the-art 
technique  and  limited  data  are  available. 

C.  JET  CHARACTERISTIC  EQUATIONS 

To  efficiently  compare  different  results,  the  widely  accepted  equations  given  below  are 
used  to  describe  round  jets  in  the  self-similar  region.  Most  modem  researchers  prefer  these 
popular  equation  formats.  Of  course,  some  conversions  were  necessary  to  adapt  all  of  the  data 
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into  this  form.  The  velocity  profiles  were  typically  approximated  using  a  Gaussian  format 
generated  from  the  jet’s  spreading  rate  as  shown  below.  The  Gaussian  profile  tends  to  be  a  good 
approximation  for  the  near  axis  regions  of  jets,  however,  the  outer  edges  of  real  jets  fall  off 
slightly  quicker  as  described  by  Malmstrom  et  al  (1997).  The  following  equations  do  not  require 
that  the  virtual  origins  to  be  co-located. 

1.  Axisymmetric  Round  Jets 


Half  Velocity  Spreading: 
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where: 


Ri/ju  -half-velocity  radius 

D  -jet  exit  diameter 

X  -axial  distance  from  the  jet  exit  nozzle 

Kiu  -velocity  spreading  rate 

Ciu  -dimensionless  velocity  geometric  virtual  origin 

Xo.iu  -velocity  geometric  virtual  origin 


Gaussian  Velocity  Profile: 
(not  all  data) 
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where: 


U  -axial  velocity 

Uc  -centerline  axial  velocity 
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(8) 


Centerline  Velocity  Decay: 
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where: 


Ue  -jet  exit  axial  velocity 

K2u  -centerline  velocity  decay  rate 

C2u  -dimensionless  velocity  kinematic  virtual  origin 

Xo,2u  -velocity  kinematic  virtual  origin 


Additional  relationships: 


R  _  T] 


(12) 


(13) 


Not  all  of  the  above  variables  were  reported  in  every  experiment.  Sometimes  the  authors 
assumed  the  virtual  origins  were  co-located  and  did  not  clearly  indicate  this  within  the  article. 
This,  in  turn,  constrains  the  reader  to  make  certain  assumptions  about  the  jet’s  characteristics, 
which  may  be  reflected  in  the  data  presented  in  Tables  1  and  2. 
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As  mentioned  before,  the  temperature  differential  between  the  surroundings  and  the  jet, 
result  in  temperature  profiles  that  are  wider  than  the  velocity  profiles.  This  same  reasoning  can 
also  be  applied  to  mass  concentration  results.  Therefore,  temperature  and  concentration  are 
interchangeable  and  characterized  in  the  following  equations 


Half  Temperature  Difference  Spreading: 


D 


(14) 


(15) 


where: 


R./,t  -half-differential  temperature  radius 
Kit  -differential  temperature  spreading  rate 
Cit  -dimensionless  temperature  geometric  virtual  origin 
Xo,u  -temperature  geometric  virtual  origin 


Temperature  Difference  Profile: 
(not  all  data) 
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where: 


T  -temperature 

Tc  -centerline  temperature 

Too  -surrounding  environment  temperature 
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Centerline  Temperature  Diff.  Decay: 


(19) 
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where: 


Ue  -jet  exit  temperature 

Kat  -centerline  differential  temperature  decay  rate 

C2t  -dimensionless  temperature  kinematic  virtual  origin 

Xo,2t  -temperature  kinematic  virtual  origin 


The  half-velocity  (temperature)  length,  Ri/2u  (Ri/2t),  represents  the  radial  distance  from 
the  jet’s  axial  centerline  to  the  position  where  the  velocity  (temperature)  is  one-half  the 
maximum  centerline  value.  The  values  of  the  dimensionless  virtual  origins  (Ciu,  C2u,  Ci„  C2t) 
and  multiplication  constants  (Kiu,  K2u,  Ku,  K2t),  in  addition  to  several  other  relevant  facts  from 
several  selected  experiments,  are  listed  in  Table  1.  As  shown,  the  multiplication  constants  (K’s) 
are  fairly  consistent  (<10%  variation)  between  experiments.  However  as  noted  by  several 
researchers,  the  virtual  origins  are  found  to  vary  widely  (-0.6  <  C’s  <  4).  Unfortunately,  no 
consistent  underlying  reason  has  yet  been  correlated  to  explain  the  wide  variations  in  the  virtual 
origin  values.  In  addition,  not  all  virtual  origins  have  been  quoted  in  the  research  which  makes 
comparisons  between  experiments  difficult. 

2.  Two-dimensional  Planar  Jets 

The  similarity  analysis  for  the  planar  jet  is,  as  expected,  very  similar  to  the  axisymmetric 
round  jet.  However,  the  characteristic  equations  take  into  account  Cartesian  (instead  of 
cylindrical)  coordinates  and  the  use  of  appropriate  substitutions. 
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The  overall  resultant  equation  forms  are: 


Half  Velocity  Spreading: 
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where:  Yi^u  -half-velocity  transverse  length 

H  -jet  height 


Gaussian  Velocity  Profile: 
(not  all  data) 
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Centerline  Velocity  Decay: 
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Additional  relationships: 


(28) 


(29) 
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The  values  of  the  dimensionless  virtual  origins  (Ci„,  C2u,  Cu,  €21)  and  multiplication 
constants  (Kiu,  K2U,  Kk,  K2t)  for  planar  jets  are  listed  in  Table  2.  Again,  the  multiplication 
constants  (K’s)  are  fairly  consistent  (<25%  variation)  between  experiments,  while  the  virtual 
origins  vary  widely  (-6H  <  C’s  <  4H). 
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D.  VARIATION  OF  EXPERIMENTAL  RESULTS 


Most  experiments  have  validated  the  similarity  analysis,  as  described  in  Chapter  H, 
Section  D.  However,  works  by  Kotsovinas  (1976)  and  Jenkins  /  Goldschmidt  (1973)  have 
continued  to  shed  doubts  on  the  linear  behavior  of  planar  jets  in  the  far  similarity  region. 
Depending  on  the  downstream  range  investigated,  the  planar  jet  centerline  decay  data  appears  to 
be  slightly  non-linear  (greater)  when  using  the  standard  similarity  equations  (26  &  35)  far  from 
the  jet  exit  (X/H  >  200).  Therefore,  the  regression  range  used  by  researchers  to  obtain  virtual 
origins  and  multiplication  constants  are  important. 

1.  Incomplete  Data 

Accurate  universal  jet  predictions  are  very  hard  to  acquire.  Several  factors  can  influence 
the  development  of  each  jet,  as  well  as  describe  their  general  shape.  Almost  no  text  or 
experimental  report  has  incorporated  all  of  the  required  factors  to  completely  characterize  jet 
flow.  Each  is  based  on  very  limited  data  and  a  general  overall  behavior  is  assumed  for  all 
comparable  jets.  Although  this  is  partially  true,  due  to  the  similarity  relationships,  one  can  not 
apply  these  results  to  all  jets  (even  when  flow  conditions  seem  to  be  the  same). 

The  most  commonly  recognized  difference  among  experiments  is  the  location  of  the 
virtual  origins.  Many  researchers  have  not  clearly  recognized  that  in  some  situations  the 
geometric  (Cm,  Cu)  and  kinematic  (C2u,  C2t)  virtual  origins  do  not  necessarily  coincide  as  would 
be  expected  by  the  similarity  arguments  [See  Goldschmidt  /  Bradshaw  (1981)  and  Flora  / 
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Goldschmidt  (1969)].  This  is  further  complicated  when  some  authors  fail  to  report  all  of  the 
virtual  origins.  This  consensus  usually  arises  from  the  assessment  that  only  the  multiplication 
constants  (K’s)  are  required  to  adequately  describe  jet  characteristics.  However,  one  can  not 
completely  describe  the  longitudinal  position  of  the  centerline  decay  and  lateral  spread  profiles 
without  knowing  the  reference  points  (virtual  origins)  upon  which  they’re  based. 

The  virtual  origin,  therefore,  is  an  important  parameter,  although  it  may  be  quite 
challenging  to  predict.  Papadopoulos  /  Pitts  (1999)  have  developed  successful  centerline 
prediction  equations  for  iso-thermal  jets  using  the  initial  momentum  flux  and  local  centerline 
turbulence  intensity.  However,  further  research  on  heated  jets  has  not  been  performed. 

2.  Jet  Variations 

To  assume  that  the  results  of  one  experiment  will  completely  match  another  is  erroneous. 
Slight  variations  in  the  measured  multiplication  constants  and  virtual  origin  values  (See  Table  1 
and  2)  can  have  significant  effects  on  jet  predictions  using  the  similarity  equations  (Equations  5- 
36). 


In  addition,  free  jets  have  been  found  to  contain  large-scale  coherent  disturbances  within 
the  shear  layers  on  top  of  the  seemingly  random  turbulent  motions  [See  Thomas  /  Prakash 
(1991),  Mourn  et  al  (1983),  Hussain  (1983),  Gutmark  /  Ho  (1983)  and  Crow  /  Champagne 
(1971)]  which  are  initial  condition  dependent.  These  disturbances  can  be  excited  relatively 
easily  and  can  drastically  effect  the  entire  motion  of  the  jet  (non-quasi-steady).  Hussein  et  al 
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(1994)  and  George  (1989)  further  described  how  a  more  general  analysis  of  the  similarity  theory 
reveals  that  a  self-preserving  state  can,  indeed,  be  uniquely  determined  by  the  initial  conditions. 
Hence,  a  family  of  farfield  similarity  solutions  could  exik  for  each  type  of  jet,  depending  on  its 

initial  conditions. 


The  density  (and/or  temperature)  dfferential  effect  on  jet  development  has  been 
investigated  by  Drobniak  et  al  (1998),  Sautet  /  Stepowski  (1995),  and  Russ  /  Stiykowski  (1993). 
Further  articles  by  Kyle  /  Sreenivasan  (1993)  and  Monkewitz  et  al  (1990)  describe  how  “hot  jet 
to  surrounding  environment”  density  ratios  less  than  approximately  pjet/p®  <  0.6  can  cause  a  self- 
excited  large-scale  coherent  mode  that  produces  a  marked  jump  in  the  spreading  rate.  This  mode 
is  caused  by  the  formation  of  a  very  regular  sequence  of  intense  vortex  ring  structures  within  the 
jet.  The  behavior  of  the  mode  has  been  found  to  be  independent  of  background  disturbances 
within  typical  experimental  Reynolds  Number  ranges  (10^2  to  10^5).  The  resulting  “strong 
pairing  process”  of  the  mode  leads  to  abnormally  large  centerline  velocity  fluctuations  and  early 
potential  core  breakup.  Current  commercial  CFD  codes,  solving  RANS  equations,  cannot  model 
these  excited  modes  which  have  become  major  obstacles  when  accurately  predicting  free  jet 

flows. 


Some  additional  experimental  variables  that  haven’t  always  been  controlled  or  reported 
include  the  following.  Fan  vibrations  felt  at  the  nozzle  and  surrounding  area  disturbances  are 
known  to  cause  significant  changes  in  the  jets,  especially  at  low  flow  rates  as  described  by 
Gutmark  /  Ho  (1983).  Modem  researchers  have  tried  to  prevent  such  interactions  to  varying 
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degrees  of  success,  but  their  reassessment  of  earlier  experiments  are  still  quite  limited.  Small  or 
semi-confmed  test  sections  and  inadequate  equipment  have  also  tainted  some  earlier  experiments 
as  described  by  Malmstrom  et  al  (1997),  Hussein  et  al  (1994),  and  Panchapakesan  (1993).  The 
overall  room  size  and  re-circulation  patterns  seem  to  deform  and  shift  the  immerging  jet  enough 
to  effect  the  results  between  different  experimental  facilities.  Extreme  care  has  been  taken  by 
modem  researchers  to  minimize  these  effects. 
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Table  1.  Axisymmetrical  Round  Jet  Experimental  Data 


25 


Table  2.  Two-dimensional  Planar  Jet  Experimental  Data 
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E.  OVERALL  TRENDS 


Several  overall  trends  in  experimental  data  have  been  noticed  throughout  the  years  with 
varying  degrees  of  consistency.  Although  a  particular  reported  trend  might  have  an  exception  in 
one  experiment,  the  overwhelming  inclinations  (as  described  by  many  additional  researchers) 
will  be  described  here. 

1.  Exit  Velocity  Effects 

Several  earlier  experiments  report  a  diameter-based  Reynolds  Number  (Red  #) 
dependence  of  jet  flows  with  exit  Rea  #’s  less  than  approximately  SxlO'*.  By  varying  the  jet  exit 
velocity,  Rea  #  dependence  was  seen  as  a  change  in  the  spreading  rate  and  centerline  decay.  A 
more  comprehensive  study  by  Malmstrom  et  al  (1997)  which  involved  varying  the  exit  velocity, 
as  well  as  the  jet  diameter,  concluded  that  the  variance  was  based  on  jet  exit  velocity,  not  the  Rea 
#.  As  jet  exit  velocity  was  increased  beyond  6  m/s,  the  velocity  centerline  decay  multiplication 
factor  (K2u)  decreased  to  an  asymptotic  0.17  value  (for  round  jets)  which  agrees  well  with  Table 
1.  Below  6  m/s,  as  jet  exit  velocity  is  decreased,  the  centerline  decay  and  spreading  rate 
increase. 


2.  Turbulence  Intensity  Effects 

The  dominant  factor  that  seems  to  effect  experimental  virtual  origin  location  is  turbulence 
intensity  at  the  jet  exit  profile.  Both  the  centerline  turbulence  intensity  and  boundary  layer 
thickness  can  drastically  shift  the  origins  in  seemingly  random  directions.  Several  researchers 
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have  attempted  to  qualify  these  effects  with  limited  success  in  narrowly  defined  jet 
configurations.  Changes  in  the  turbulence  intensity  have  also  sometimes  masked  the  effects  of 
velocity  dependence  on  jet  development.  As  the  jet  exit  velocity  increases,  the  core  turbulence 
intensity  generally  drops  (unless  artificially  held  constant)  as  the  boundary  layers  become 
thinner.  Therefore,  as  exit  turbulence  intensity  decreases,  the  jet  spreads  less  and  the  centerline 
decay  decreases.  However,  Goldschmidt  /  Bradshaw  (1981),  Kleis  /  Floss  (1974),  and  Flora  / 
Goldschmidt  (1969)  have  shown  that  the  above  Tu  effects  still  hold  true  when  the  exit  velocity  is 
held  constant. 

F.  TARGET  VALUES 

In  order  to  compare  the  CFD  results  to  experiments,  target  values  must  be  chosen.  A 
band  (or  range)  of  experimental  data  was  select  over  a  specific  experiment  because  every  jet 
tends  to  be  unique  as  described  in  Section  (D).  After  analyzing  Tables  1  and  2,  evaluating  the 
references,  and  studying  more  recent  investigations,  the  following  parameters  were  chosen  to 
compare  and  model  typical  jet  flow  for  this  research. 

1.  Centerline  Decay  and  Spreading  Values 

As  shown  in  Tables  1  and  2,  the  spreading  rates  and  centerline  decay  values  are  fairly 
constant,  but  do  exhibit  a  slight  variation  (-10%  for  the  round  jet  and  -25%  for  the  planar  jet). 
Based  on  measurement  technique,  data  scatter,  test  range,  researcher  credibility,  the  following 
values  will  be  used  to  compare  CFD  results. 
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Table  3.  CFD  Simulation  Target  Values 


Flow  Type 

Vel 

Spread 

Rate 

Geom. 

Virt. 

Orig. 

C/L 

Vel. 

Decay 

Kin. 

Virt 

Orig. 

Temp 

Spread 

Rate 

Spread 

Virt. 

Orig. 

C/L 

Temp 

Decay 

Decay 

Virt 

Orig. 

|K,„) 

{C,„} 

IKzul 

{C2„} 

(Kit) 

(Cit) 

{K2t} 

{C2t} 

Round  Jet 

0.085- 

0.095 

-1.0  to 
4.0 

0.160- 

0.172 

-1.0  to 
4.0 

0.098- 

0.110 

0.0  to 
2.0 

0.185- 

0.194 

0.0  to 
2.0 

Var. 

11.7% 

5D 

7.5% 

5D 

12.2% 

2D 

4.9% 

2D 

Planar  Jet 

0.096- 

0.120 

-6.0  to 
2.0 

0.155- 

0.220 

-2.0  to 
4.0 

0.125- 

0.140 

-5.0  to 
0.0 

0.250- 

0.280 

-5.0  to 
0.0 

Var. 

25% 

8H 

41.9% 

6H 

12% 

5H 

12% 

5H 

Target  virtual  origin  ranges  were  selected  to  include  the  extreme  reported  variations.  As 
mentioned,  current  steady-state  Favre-Averaged  Navier-Stokes  (FANS)  CFD  codes  can  not 
model  many  important  external  factors  (i.e.,  nozzle  vibration,  room  disturbances,  etc.)  that  can 
significantly  effect  the  virtual  origin  position.  Therefore,  the  virtual  origin  target  range  is  used 
only  for  vague  simulation  result  comparisons. 

2.  Turbulent  Prandtl  Number 

The  turbulent  Prandtl  Number  (Prt)  quantifies  the  transport  rate  of  momentum  to  that  of 
heat  and  was  used  as  an  input  to  the  CFD  code  k-e  turbulence  /  energy  model.  Free  jet  flows, 
which  are  not  effected  by  walls,  display  different  Prt  variations  normal  to  the  primary  flow 
direction  than  typical  internal  flows  (i.e.  equilibrium  pipe  flow).  The  turbulent  Prandtl  Number 
varies  across  the  round  jet’s  profile  and  was  quantified  by  Reynolds  (1976)  as  shown  in  Table  4 
and  Figure  4. 
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Table  4.  Turbulent  Prandtl  Numbers  for  Jets 


Flow  Type 

Local  Prt 

Recommended 
Global  Prt 

ys=l 

7^0.5 

y-^0 

Round  Jet 

0.73 

1.2 

1.7 

0.7 

Planar  Jet 

0.5 

0.9 

1.3 

0.5 

Figure  4.  Typical  Prt  and  y  Variations  across  a  round  and  planar  jets 

[Based  on  Reynolds  (1976)] 
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The  Turbulence  intermittency  factor  (y)  in  the  jet  profile  varies  from  unity  near  the  jet 
centerline  (<L)  and  inner  part  of  its  shear  layer  to  0.5  and  below  in  the  outer  edges  of  the  shear 

layer.  The  most  striking  feature  of  the  results  is  the  difference  between  the  Prt  in  the  body  (yi=l) 
and  the  outer  edges  of  the  highly  turbulent  shear  layer.  A  Pr,  greater  than  unity  indicates  that 
momentum  is  being  transferred  more  readily  than  heat  in  these  highly  intermittent  regions  via 
relatively  large-scale  turbulent  structures. 

These  variations,  although  significant  across  the  whole  jet,  are  fairly  small  across  the 
largest  inner  portion  of  the  jet  and  can  be  approximated  using  a  single  global  value  for  the  entire 
jet.  In  general,  the  body  values  for  round  jets  are  larger  than  for  planar  jets.  Chen  /  Rodi  (1980) 
and  Reynolds  (1976)  recommend  the  well  accepted  global  Pr,‘s  listed  in  Table  4  when  the 
density  differences  are  small  and  the  molecular  Prandtl  Number  (Pr)  is  close  to  unity  (Prt, air  ~ 
0.7).  These  global  approximations,  of  coarse,  result  in  small  CFD  temperature  over-prediction 
errors  in  the  outer  edges  of  the  jets. 
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IV.  CFD  SIMULATION 


A.  CODE  AND  COMPUTER 

The  CFD  Research  Corporation  (CFDRC)  developed  the  computational  fluid  dynamics 
package  used  in  this  research  to  analyze  the  different  jet  configurations.  The  CFD-ACE+ 
program  (version  6.2)  was  run  on  an  EBM/PC-based  Gateway  (Model  GP6-450)  450  MHz  clock- 
speed  computer.  The  Random  Access  Memory  was  upgraded  to  384MB  to  handle  grid  files  with 
up  to  300,000  three-dimensional  unstructured  cells.  A  brief  description  of  the  CFDRC  code  and 
its  operation  is  provided  in  Appendix  A. 

B.  THE  PROCESS 

To  gain  experience  with  the  code,  isothermal  jets  were  initially  modeled  using  small 
domains  with  a  coarse  grid  to  minimize  solution  run  times.  This  process  allowed  the  flow 
conditions  to  be  fully  understood  before  the  complexity  of  heat  transfer  was  added.  After  the 
code  was  mastered  at  each  step,  an  additional  feature  ,  or  component  was  activated.  Grid 
refinement  was  also  developed  in  the  high  gradient  areas.  Only  the  final  heated  jet  models  are 
discussed  within  this  report.  In  the  anticipation  of  further  research  with  more  complex  flows 
using  this  code,  only  3-dimensional  unstructured  tetrahedral  cells  were  used  in  the  generation  of 
all  grid  domains. 
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c. 


HEATED  AXISYMMETRICAL  JET  MODEL 


The  heated  axi symmetrical  round  jet  was  modeled  to  investigate  the  ability  of  the  code  to 
predict  these  types  of  flows.  Although  this  code  can  not  model  all  factors  effecting  jet 
development  (i.e.,  coherent  structural  modes,  etc.),  an  overall  performance  evaluation  of  the  code 
compared  to  experimental  data  was  desired. 

1.  Grid  Domain 

A  one-quarter  section  axisymmetrical  grid  was  used  to  accurately  model  the  developing 
jet  and  its  surroundings  while  minimizing  the  total  grid  size.  Smaller  sections  (i.e.  1/6,  1/8,  etc.) 
were  not  used  to  prevent  possible  symmetric  edge  interference  problems.  Typical  grid 
dimensions  and  overall  configuration  are  shown  in  Figure  5.  The  total  grid  length  was  selected 
to  be  40  jet  exit  diameters  with  a  grid  radius  of  slightly  less  than  8  jet  exit  diameters  in  order  to 
capture  the  flow  features  without  boundary  interference.  The  overall  grid  domain  was  divided 
into  two  separate  sub-domains.  The  inner  sub-domain  consisted  of  smaller  cells  to  more 
accurately  capture  the  high  gradients  of  the  main  jet  flow  and  its  free  shear  layer  (adaptive 
meshing  was  not  available).  The  inner  sub-domain  was  constmcted  as  an  inverted  truncated 
cone  to  mimic  the  jet’s  radial  spread  as  it  travels  downstream.  The  inner  sub-domain 
encompasses  most  of  the  jet’s  flow  throughout  the  entire  grid  domain.  The  jet  enters  this  sub- 
domain  through  the  lower  face  and  exits  through  the  top  with  entrainment  coming  in  from  the 
sides.  The  bottom  boundary  of  the  inner  sub-domain  is  slightly  larger  than  the  jet  exit  diameter 
to  facilitate  a  finer  mesh  (more  accurate  solution)  around  this  transitional  flow  area. 
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Figure  5.  AsSsymmetrica!  Grid 


Before  generating  the  internal  tetrahedral  volume  grid,  the  inner  sub-domain’s  botmdary 
surface  grids  were  automatically  generated  using  the  following  dimensionless  control 
parameters. 

•  Maximum  cell  size:  0.007 

•  Minimum  cell  size:  1 . 143x  1 0'^ 

•  Curvamre  Criterion:  8 

•  Surface  Transition  Factor:  1 . 1 
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The  Curvature  Criterion  defines  the  approximate  maximum  angle  of  curvature  per  cell 
that  will  be  generated  at  edges  that  are  curved.  The  Transition  Factor  limits  the  rate  of  cell 
growth  (size  change)  when  generating  triangular  surface  cells  from  a  curved  edge  to  areas  far 
from  all  curved  edges. 

After  all  surfaces  of  the  inner  sub-domain  were  gridded,  a  separate  volume  tool  was  used 
to  invoke  the  tetrahedral  volume  grid  generation  algorithm.  The  tetrahedral  volume  grid 
generation  algorithm  used  in  CFD-GEOM  is  a  variation  of  the  advancing  front  grid  generation 
technique.  The  size  of  the  generated  tetrahedrons  change  smoothly,  based  on  the  properties  of 
the  surface  grids.  One  additional  mechanism  (besides  those  used  in  surface  generation)  was  used 
to  control  the  tetrahedral  grid  volume  generation.  The  Volume  Transition  Factor,  similar  to  the 
Surface  Transition  Factor,  determined  how  far  off  a  face  that  a  field  point  would  be  generated  for 
the  formation  of  a  new  tetrahedron.  For  the  inner  sub-domain,  the  Volume  Transition  Factor  was 
set  at  1.1 

In  general,  the  outer  sub-domain  contained  larger  cells  to  model  the  airflow  in  the 
relatively  slower  moving  surrounding  environment.  The  outer  boundaries  of  this  domain  were 
made  vertical  on  the  sides  and  horizontal  on  the  top  and  bottom  to  simplify  boundary  conditions 
at  these  areas.  The  grid  generation  controls  for  this  domain  were  set  as  follows. 
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•  Maximum  cell  size:  0.05 


•  Minimum  cell  size:  1.143x10'^ 

•  Curvature  Criterion:  30 

•  Surface  Transition  Factor:  1.2 

•  Volume  Transition  Factor:  1.2 

The  outer  sub-domain  grid  was  generated  after  the  inner  sub-domain  so  that  the 
tetrahedral  cells  could  progress  smoothly  from  the  small  inner  cells  to  the  larger  outer  sub- 
domain  cells  at  the  grid  boundaries. 

Overall,  the  entire  domain  contained  317,273  tetrahedral  cells  with  the  following  grid 
qualities  as  computed  by  the  code.  The  CFDRC  User  Manual  (2000)  can  be  consulted  for 
further  definition  of  each  criterion. 

•  Dihedral  Angle:  3  cells  were  less  than  5  degrees  (default  min.=5  degrees) 

•  Worst  Centroid/Face  Angle:  7.24  degrees  (default  min.=5  degrees) 

•  Sliver  Quality:  5  cells  were  greater  than  7  (default  max.=7) 

•  Skew  Quality:  52  cells  were  greater  than  0.95  (default  max.=0.95) 

Only  a  small  faction  of  the  cells  in  non-critical  areas  violated  the  default  limits  and  the 
gridded  domain  was  determined  to  be  satisfactory.  Several  other  complex  domain  configurations 
were  tried  to  further  concentrate  small  cells  in  high  gradient  areas  without  greater  grid  quality 
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success.  In  addition,  the  higher  domain  complexity  adversely  effected  the  solution  convergence. 
More  controls  are  needed  to  properly  manage  the  automatic  grid  generation  process  for 
unstructured  cells  and  CFDRC  plans  to  release  future  versions  addressing  this  problem. 

2.  Solver  Setup 

a.  Overall 

Table  5  below  provides  a  summary  of  all  solver  [CFD-GUI  /  CFD-ACE  (U)] 
inputs.  To  simulate  the  heated  jet,  the  flow,  heat  transfer  and  turbulence  modules  were  activated. 
The  relatively  simple  k-z  turbulence  model  was  chosen  to  gain  experience  manipulating  the  code 
while  maximizing  solution  convergence.  As  mentioned  in  Appendix  B,  the  k-t  turbulence  model 
constants  can  not  be  changed  (within  the  code)  for  axis3mimetrical  flows  as  suggested  by 
Launder/Spalding  (1974).  Therefore,  without  being  able  to  vary  and  Cei  (turbulence  model 
closure  constants)  across  the  profile,  the  overall  prediction  accuracy  for  axisymmetrical  jets  was 
expected  to  be  poor. 

b.  Volume  Conditions 

Air  at  ambient  conditions  was  selected  as  the  computational  fluid.  To  simulate 
the  heated  jet  (AT  =  30K),  density  calculations  evaluated  using  the  Ideal  Gas  Law  while  the 
viscosity  was  determined  by  Sutherland’s  Law  [|i=AT^^^/(B+T)],  where  A  and  B  are  constants 
specified  in  Table  5.  Since  the  simulation  involved  small  differential  temperatures,  fluid 
conductivity  {K)  and  specific  heat  (Cp)  evaluations  were  determined  by  selecting  the  constant 
molecular  Prandtl  Number  (Pr  =  0.7)  option  and  constant  Cp  =  1000  J/kg-K,  respectively. 
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c.  Boundary  Conditions 

The  grid  inlet  boundary  was  modeled  as  0.0254m  (1-inch)  diameter  jet  exit  nozzle 
with  a  constant  velocity  (top  hat)  profile.  Since  the  grid  domain  was  modeled  as  a  1/4*  section, 
the  inlet  boundary  is  shaped  like  a  1/4*  “pie  section”  with  a  radius  D/2  =  0.0127m.  The  ideal 
‘top  hat’  outlet  jet  velocity  profile  approximates  the  flow  that  most  researches  try  to  develop 
using  large  nozzle  contractions.  Although  ‘real’  jets  have  small  boundary  layers  which  can 
drastically  effect  their  development,  an  ideal  top  hat  profile  was  used  to  simplify  the  CFD  model. 
With  a  properly  supplied  sub-routine,  a  user  can  define  any  desirable  inlet  jet  velocity  profiles, 
however,  the  improvement  in  the  final  CFD  solution  was  expected  to  be  negligible  when 
compared  to  the  widely  varying  experimental  jet  data. 

The  constant  jet  nozzle  exit  velocity  (Ue)  {grid  inlet}  was  initially  set  at  40m/s,  in 
the  axial  direction  (x)  ,  with  a  uniform  temperature  of  330  degrees,  K.  All  inlet  turbulence 
kinetic  energies  (k)  needed  for  the  Favre-Averaged  Navier-Stokes  (FANS)  conservation 
equations  (Equations  42  and  43)  were  calculated  from  Equation  37,  in  which  the  turbulence 
intensity  (Tm)  was  assumed  to  be  contemporary  experimental  values.  A  turbulence  intensity  of  Tu 
=  0.5%  was  originally  selected  which  resulted  in  A:  =  0.06  mW. 

(37) 
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The  inlet  dissipation  rate  (s)  was  input  into  the  code  by  using  the  ‘length  scale 
option’  as  discussed  in  Appendix  B.  This  method  required  that  the  inlet  size  dimension  (L)  be 
specified,  then  the  turbulent  length  scale  (0.3%  of  the  size  of  the  inlet  size  dimension)  is 
subsequently  calculated  within  the  program.  Since  the  grid  inlet  boundary  is  a  “pie  section”,  the 
size  dimension  was  set  at  L  =  0.0127m  (inlet  boundary  radius). 

The  straight  sides  of  the  1/4*  section  grid  domain  were  set  as  symmetric 
boundaries.  The  rest  of  the  grid  boundaries  (top,  curved  front  surface  and  rest  of  the  bottom) 
were  established  as  fixed  (static)  ambient  pressure  outlets.  This  setup,  however,  does  allow  flow 
to  enter  the  grid  domain  through  the  side  and  bottom  boundaries  for  entrainment.  The  outlet 
boundary  parameters  were  set  to  characterize  the  flow  as  it  entered  the  grid  domain  from  the 
infinite  surroundings. 

d.  Initial  Conditions 

The  initial  conditions  are  displayed  in  Table  5.  The  starting  turbulent  viscosity 
was  artificially  elevated  (via  k  and  L)  to  enhance  solution  convergence  during  the  earliest 
solution  iterations. 


e.  Solver  Control 

A  l®'-order  upwind  spatial  differencing  scheme  was  chosen  for  all  variables  in  the 
jet  models  to  simplify  convergence  operations  and  minimize  solution  times  (typically  2-3  days). 
The  default  CGS-nPre  (conjugate-gradient-squared  plus  preconditioning)  linear  equation  solver 
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was  also  chosen  to  resolve  the  algebraic  equations  for  each  dependent  variable  and  the 
controlling  parameters.  The  solver  performed  satisfactory  without  approaching  its  internal 
default  maximum  number  of  sweeps  criterion  (50  or  500  sweeps,  depending  on  the  variable). 
The  under-relaxation  parameters  were  also  left  at  their  default  values  (see  Table  5),  which 
seemed  to  work  adequately. 
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Modules 


Outlets 


Heat 


Turbulence 


Table  5.  Summary  of  Round  Jet  Solver  Settings 


_  Problem  Type 


Flow 

Heat  Transfer 

_  T  urbulence 


_  Steady  Time  Dependence 


_  Model  Options 


Flow 

Reference  Pressure 

Turbulence 

K-Epsilon  Model 

Pn  =  0.7 


Volume  Conditions 


Both 

Property  Mode 

Fluid 

Domains 

Density  Evaluation  Method 

Ideal  Gas  Law 

M.W.  =  29  (default) 

Viscosity  Evaluation  Method 

Sutherland’s  Law 

A  =  1.4605E-6 

B  =  112  (default) 

Specific  Heat  Evaluation  Method 

Constant 

Cp  =  1 000  J/ka-K  (default) 

Conductivity  Evaluation  Method 

Mol.  Prandtl  Number 

Pr  =  0.7 

Boundary  Conditions 


Mode 
U  Veloci 
V  Veloci 
W  Veloci 

Relative  pressure _ 

Temperature _ 

Turbulent  Kinetic  Ener 
Dissipation  Rate  fLenqth  Scale  Method 


Mode _ 

Relative  pressure _ 

Temperature _ 

Turbulent  Kinetic  Ener 
Dissipation  Rate  f Length  Scale  Method 


_  Initial  Conditions 


Initial  Source 
U  Veloci 
VVeloci 
W  Veloci 
Relative  pressure 


Temperature 


Turbulence  Kinetic  Ener 
Dissipation  Rate  (Length  Scale  Method 


Solver  Control 


Fix  VeL  (Cartesian 


40  (m/s 


Spatial  Differencing  Scheme 


Solver  Selection 


Under-Relaxation  Parameters 


I  Velocities 

I  Enthal 
Turbulence 


Veloci 
P  Correction 
Enthalpy 
Turbulence 


Velocities 
P  Correction 
Enthalpy 
Turbulence 
Pressure 

Viscosity 

Temperature 


0.0127  (m 


Fixed  Pressure 


Constant:  0.00015  (m  /s 


Constant  Values 


10  (m/s 


0  (m/s 


0.01  (m 


Upwind  (default 


Upwind  (default 


Upwind  (default 


Upwind  (default 


CGS+Pre  (default 


CGS+Pre  (default 


CGS+Pre  (default 


CGS+Pre  (default 


0,2  (default 


0.2  (default 


0.05  (default 


0.2  (default 


1  (default 


1  (default 


1  (default 
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D.  HEATED  PLANAR  JET  MODEL 


1.  Grid  Domain 

Again,  only  3-dimensional  unstructured  tetrahedral  cells  were  used  in  the  generation  of 
the  planar  grid  domain.  A  rectangular  grid  was  used  to  model  half  (1/2)  of  the  developing  planar 
jet  and  its  still  surroundings.  See  Figure  6  below  for  typical  grid  dimensions  and  overall 
configuration.  The  total  grid  length  was  selected  to  be  40  jet  exit  heights  (H).  A  grid  width  of 
approximate  11.5  jet  exit  heights  and  depth  of  4  jet  exit  heights  was  established  in  order  to 
capture  the  flow  features  without  boundary  interference.  The  overall  grid  domain  was  again 
divided  into  two  separate  sub-domains.  The  inner  sub-domain  consisted  of  smaller  cells  to  more 
accurately  capture  the  high  gradients  of  the  main  jet  flow  and  its  free  shear  layer.  The  jet  enters 
the  inner  sub-domain  through  the  lower  face  and  exits  through  the  top  with  entrainment  coming 
from  the  non-symmetrical  far  side.  The  bottom  boundary  of  the  inner  sub-domain  is  slightly 
larger  than  the  jet  exit  diameter  to  facilitate  a  finer  mesh  (more  accurate  solution)  around  this 
transitional  flow  area. 
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Figure  6.  Two-Dimensional  Jet  Grid 

The  inner  sub-domain’s  boundary  surface  and  volume  grid  was  automatically  generated 
using  the  following  control  parameters. 

•  Maximum  cell  size:  0.004 

•  Minimum  cell  size:  5.2x10'^ 

•  Surface  Transition  Factor:  1 . 1 

•  V  olume  T  ransition  Factor:  1 . 1 
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The  outer  sub-domain  grid  generation  controls  were  as  follows. 


•  Maximum  cell  size:  0.01 

•  Minimum  cell  size:  5.2x10’® 

•  Surface  Transition  Factor:  1.1 

•  Volume  Transition  Factor:  1 . 1 

The  outer  sub-domain  grid  was  generated  after  the  inner  sub-domain  so  that  the 
tetrahedral  cells  could  progress  smoothly  from  the  small  inner  cells  to  the  larger  outer  sub- 
domain  cells  at  the  grid  boundaries.  The  typical  cell  size  was  smaller  than  the  axisymmetric  jet 
model  due  to  the  smaller  overall  grid  volume  (limited  to  approximately  300,000  cells). 

Overall,  the  entire  domain  contained  297,657  tetrahedral  cells  with  the  below  grid 
qualities  as  computed  by  the  code.  All  grid  quality  checks  were  within  limits  and  the  gridded 
domain  was  determined  to  be  satisfactory. 

•  Worst  Dihedral  Angle:  6.54  degrees  (default  min.=5  degrees) 

•  Worst  Centroid/Face  Angle:  18.1 1  degrees  (default  min.=5  degrees) 

•  Sliver  Quality:  1  cell  at  7.01  (default  max.=7) 

•  Worse  Skew  Quality:  0.93  (default  max.=0.95) 


45 


2. 


Solver  Setup 


a.  Overall 

Table  6  below  provides  a  summary  all  of  the  solver  inputs,  which  are  very  similar 
to  the  round  jet  simulation. 

b.  Boundary  Conditions 

The  2-dimensional  jet  exit  nozzle  (grid  inlet  boundary)  was  modeled  with  a  height 
of  0.013m  (1/2-inch)  and  a  constant  ‘top  hat’  velocity  profile.  Since  the  grid  domain  was  a  1/2 
section,  the  inlet  boundary  has  a  half-height  of  D/2  =  0.0065m.  Therefore,  the  inlet  dissipation 
rate  (e),  using  the  length  scale  option,  was  set  at  L  =  0.0065m  to  match  the  grid  inlet  size  (similar 
to  the  axis)mimetrical  jet  model). 

The  bottom  boundary  was  established  as  an  adiabatic  wall  to  emulate  the 
configuration  of  most  heated  planar  jet  experiments.  The  rest  of  the  grid  boundaries  were 
established  as  fixed  ambient  pressure  outlets.  The  outlet  boundary  parameters  were  again  set  to 
characterize  the  surrounding  flow  if  it  entered  the  grid  domain  from  the  infinite  surroundings. 

c.  Initial  Conditions 

The  initial  conditions  are  displayed  in  Table  6.  The  starting  turbulent  viscosity 
was  artificially  elevated  (via  k  and  L)  to  enhance  solution  convergence  during  the  earliest 
iterations. 


46 


d.  Solver  Control 


A  l®‘-order  upwind  spatial  differencing  scheme  was  chosen  for  all  variables  in  the 
jet  models  to  simplify  convergence  operations  and  minimize  solution  times  (typically  3  days). 
The  default  CGS+Pre  (conjugate-gradient-squared  plus  preconditioning)  linear  equation  solver 
was  also  chosen  to  resolve  the  algebraic  equations  for  each  dependent  variable  and  the 
controlling  parameters.  The  solver  performed  satisfactory  without  approaching  its  internal 
default  maximum  number  of  sweeps  criterion  (50  or  500  sweeps,  depending  on  the  variable). 
Again,  the  under-relaxation  parameters  were  also  left  at  their  default  values,  which  seemed  to 
work  adequately. 


47 


Table  6. 


Summary  of  Planar  Jet  Solver  Settings 


Modules 

Flow 

Heat  Transfer 

Turbulence 

Global 

Steady  Time  Dependence 

Model  Options  | 

Flow 


Turbulence 


Heat 


Turbulence 


Reference  Pressure 


K*Epsilon  Model 


Pr,  =  0.5 


Volume  Conditions 


Both 

Fluid 

Domains 

Density  Evaluation  Method 

Ideal  Gas  Law 

M.W.  =  29  (default) 

Viscosity  Evaluation  Method 

Sutherland’s  Law 

A  =  1.4605E-6 

B  =  112  (default) 

Specific  Heat  Evaluation  Method 

Constant 

Cp  =  1000J/kg-K  (default) 

Conductivity  Evaluation  Method 

Mol.  Prandtl  Number 

Pr  =  0.7 

EBl 


Boundary  Conditions 


Mode 


U  Veloci 


V  Veloci 


W  Veloci 


Relative  pressure 


Temperature 


Turbulent  Kinetic  Ener 


Dissipation  Rate  {Lenqth  Scale  Method 


Mode 


Relative  pressure 


Temperature 


Turbulent  Kinetic  Ener 


Dissipation  Rate  {Length  Scale  Method 


Initial  Conditions 


Initial  Source 


Fix  Vel.  {Cartesian 


IIUS91 


1^ 


V  Veloci 


W  Veloci 


Relative  pressure 


Temperature 


Turbulence  Kinetic  Enen 


Dissipation  Rate  {Length  Scale  Method 


Solver  Control 


300  {K 


Constant  0.00015  (m  /s 


0.05 


Constant  Values 


10  (m/s 


0  (m/s 


Spatial  Differencing  Scheme 


Solver  Selection 


i  Velocities 

I  Enthal 
Turbulence 


P  Correction 

Enthalpy 

Turbulence 


Upwind  (default 


Upwind  (default 


Upwind  (default 


CGS+Pre  (default 


Temperature 


V.  COMPARISON  OF  CFD  RESULTS  TO  EXPERIMENTS 


CFD  simulation  velocity  and  temperature  profiles  at  different  axial  distances  for  both  jets 
are  displayed  in  Figures  7-8  for  the  round  jet  and  Figures  15-16  for  the  planar  jet.  By  using  the 
similarity  equations  (Equations  5  through  36),  the  non-dimensional  plots  (Figures  9-14  for  round 
jet  and  Figures  17-22  for  planar  jet)  were  generated  to  investigate  the  self-similar  ranges  of  each 
jet.  As  shown,  the  simulated  round  jet  doesn’t  achieve  self-similarity  until  approximately  30 
diameters  downstream  of  the  jet  exit  (see  Figures  9, 10,  13,  and  14).  The  simulated  planar  jet,  on 
the  other  hand,  reached  similarity  within  approximately  25  nozzle  heights  (see  Figures  17,  18, 
21,  and  22). 

The  centerline  decay  and  spreading  rates,  as  well  as  the  virtual  origins,  for  planar  and 
round  jet  simulations  are  listed  in  Table  7.  Regression  ranges  were  calculated  in  the  similarity 
region  between  30  and  40  jet  diameters  (or  slot  heights)  for  both  jet  types.  Several  values  (half¬ 
radius/length  or  centerline)  were  analyzed  by  a  Microsoft  Excel  spreadsheet  regression  program 
to  determine  the  slope  [spreading  or  centerline  decay  rate  (K)]  and  imaginary  X-intercept  [virtual 
origin  (C)]  of  each  plot.  See  Appendix  C  for  a  typical  sample  of  the  Excel  spreadsheet  used  for 
round  jet  calculations. 
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o  X=0.127m  ()CO=5) 

A  X=0.254m  (XO=10) 
X  X=0.508m  (X/D=20) ' 
o  X=0.762m  (X/D=30) 
+  )«fc0.889m  (X/D=35) , 

- X=1.016m  (X/D=40) 

o  )fe:1.143m  (X/D=45) 
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figure  8.  Round  Jet  Temperature  Profiles 


Figure  10.  Round  Jet  Centerline  Temperature  Decay 
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Figure  11.  Round  Jet  Velocity  Profile  in  the  Self-similar  Region 


0  0.05  0.1  0.15  0.2  0.25  0.3 


Eta=R/(X-C1tD) 


Figure  12.  Round  Jet  Temperature  Profile  in  the  Self-similar  Region 
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Eta=Y/(X-C1uH) 

Figure  13.  Round  Jet  Velocity  Similarity 


Eta=Y/(X-C1tH) 

Figure  14.  Round  Jet  Temperature  Similarity 
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Y(m) 

Figure  15.  Planar  Jet  Veloci^  Profiles 


Y(m) 


Figure  16.  Planar  Jet  Temperature  Profiles 


Figure  17.  Planar  Jet  Centerline  Velocity  Decay 


Figure  18.  Planar  Jet  Centerline  Temperature  Decay 
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Figure  19.  Planar  Jet  Velocity  Profile  in  the  Self-similar  Region 


Figure  20.  Planar  Jet  Temperature  Profile  in  ttie  Self-similar  Region 
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(T-TinO/(Tc-Tlnf)  U/Uc 


Figure  21.  Planar  Jet  Velocity  Simflarity 


Figure  22.  Planar  Jet  Temperature  Similarity 
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Table  7.  CFD  Simulation  Predictions 


Flow  Type 

Vel. 

Spread 

Rate 

Geom. 

Virt. 

Orig. 

C/L 

Vel. 

Decay 

Kin. 

Virt. 

Orig. 

Temp. 

Spread 

Rate 

Spread 

Virt. 

Orig. 

C/L 

Temp. 

Decay 

Decay 

Virt. 

Orig. 

{K,„} 

{C,u} 

{K2J 

{C2J 

{Ki.} 

{Ch} 

{K2.} 

{C2tl 

Round 

CFD 

0.112 

9.3 

10.0 

0.134 

11.1 

0.247 

12.3 

Jet 

Target 

Values 

0.085- 

0.095 

10^9 

Eval. 

+24% 

+7.8D 

+21% 

+7.5D 

+29% 

+10.1D 

+30% 

+11.3D 

Planar 

CFD 

0.092 

3.7 

5.2 

0.149 

$.6 

0.209 

9.1 

Jet 

0.096- 

0.120 

-6.0  to 
2.0 

0.125- 

0.140 

-5.0  to 
0.0 

Eval 

-17% 

+5.7H 

-22% 

+6.2H 

+12% 

+11.1H 

-27% 

+11.6H 

The  CFD  predictions  can  be  compared  to  the  target  values  listed  in  Table  7  and 
visualized  in  Figures  9  through  12  for  the  round  jet  and  17  though  20  for  the  planar  jet.  These 
non-dimensionalized  CFD  profile  plots  were  generated  at  X/D=40  for  the  round  jet  (X/H=35  for 
the  planar  jet).  The  target  value  ranges,  indicated  on  the  figures  for  each  jet,  represent  the 
reported  variation  of  experimental  data  (target  value  range)  in  the  self-similar  region  only.  The 
‘High  Target’  curves  indicate  the  upper  bound  on  experimental  data  while  the  ‘Low  Target’ 
curves  represent  the  lower  experimental  data  boundaries.  As  indicated,  none  of  the  CFD 
predictions  fell  within  the  targeted  ranges.  Evaluations  were  performed  to  quantify  the 
difference  between  CFD  results  and  the  target  value  ranges.  Each  CFD  value  (K  or  C)  was 
compared  to  the  ‘middle  value’  of  the  target  value  range  and  presented  as  the  ‘Eval’  values  in 
Table  7. 
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A.  HEATED  AXISYMMETRICAL  JET  RESULTS 


As  expected,  the  jt-£  turbulence  model  simulated  the  round  jet  poorly  as  discussed  in 
Section  IV.C.2  and  Appendix  B.  The  round  jet  results  indicated  that  the  code  over-predicted  the 
velocity-spreading  rate  by  24%  and  the  temperature-spreading  rate  by  29%  when  compared  to 
the  middle  of  their  respective  target  value  range.  In  addition,  the  centerline  velocity  and 
temperature  decay  rates  were  also  over-predicted  by  21%  and  30%,  respectively.  The  geometric 
and  kinematic  virtual  origins  were  over-predicted,  as  well,  by  approximately  7.5  diameters  for 
the  velocity  profiles  and  10.5  diameters  for  the  temperature  profiles.  Although  experimental 
virtual  origin  data  display  significant  variation,  typical  values  are  usually  reported  close  to  the  jet 
exit  [-1  <  X/D  <  4]. 

B.  HEATED  TWO-DIMENSIONAL  JET  RESULTS 

Planar  jet  simulation  was  generally  closer  to  experimental  data  ranges,  with  an  under¬ 
prediction  of  the  velocity-spreading  rate  of  approximately  17%  with  an  over-predicted 
temperature-spreading  rate  of  12%.  The  centerline  velocity  and  temperature  decay  rates  were 
both  under-predicted  at  22%  and  27%,  respectively.  Again,  the  geometric  and  kinematic  virtual 
origins  were  over-predicted  by  approximately  7.5  slot  heights  for  the  velocity  profiles  and  10.5 
slot  heights  for  the  temperature  profiles. 
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C.  DIMENSIONAL  JET  COMPARISONS 


To  evaluated  the  actual  differences  between  the  CFD  results  and  experiments,  the  non- 
dimensional  experimental  parameters  (K’s  and  C’s)  were  used  to  generate  actual  profiles 
[velocity  (m/s)  and  temperature  (K:  degrees  Kelvin)]  using  the  similarity  equations  (Equations  5 
through  36).  By  comparing  the  actual  profiles,  further  insight  was  discovered  about  the  accuracy 
of  the  CFD  code. 

The  ‘middle  values’  of  the  target  value  ranges  (experimental  data  ranges)  were  chosen  for 
comparison  with  the  CFD  profiles  (see  Table  8).  These  middle  values  were  used  to  generate 
‘average  experimental  profiles’  and  were  not  based  on  any  specific  study,  but  did  establish  a 
foundation  for  evaluation. 


Table  8.  Simulation  Comparison  Parameters 


Flow  Type 

Vel. 

Spread 

Rate 

Geom. 

Virt. 

Orig. 

C/L 

Vel. 

Decay 

Kin. 

Virt. 

Orig. 

Temp. 

Spread 

Rate 

Spread 

Virt. 

Orig. 

C/L 

Temp. 

Decay 

Decay 

Virt 

Orig. 

{Ki„} 

{Clu} 

|K2u} 

(C^J 

{K,t} 

{C,.} 

{K2.} 

{C,.} 

Round 

CFD 

0.112 

9.3 

0.201 

10.0 

0.134 

11.1 

12.3 

Jet 

Middle 

Target 

0.090 

■B 

0.166 

m 

Planar 

CFD 

0.092 

3.7 

5.2 

0.149 

8.6 

0.209 

9.1 

Jet 

Middle 

Target 

0.108 

0.188 

1.0 

0.133 

-2,5 

.265 

-2.5 
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The  CFD  jet  simulation  exit  velocities  [Ue=40m/s  round  jet  (Ue=35m/s  planar  jet)], 
temperatures  (Te=330K),  and  jet  sizes  [D=0.254m  round  jet  (H=0.013  planar  jet)]  were  used  as 
inputs  to  the  produce  the  average  experimental  profiles.  This  approach,  of  coarse,  ignores  the 
effects  of  velocity  and  temperature  on  the  development  of  coherent  stractures  within  free  jets. 
However,  an  overall  comparison  was  still  desired. 


The  CFD  and  average  experimental  profiles  were  evaluated  at  X/D  =  40  for  the  round  jet 
and  are  displayed  in  Figures  23  and  24.  The  higher  CFD  centerline  values  (Uc,  Tc)  are  the  result 
of  lower  Ue/Uc  and  ATe/ATc  values  at  X/D  =  40  as  displayed  in  Figures  9  and  10.  Even  though, 
the  CFD  results  displayed  higher  centerline  decay  rates  (K2u,  Kat),  their  excessively  large  virtual 
origins  (C2u,  Cat)  produced  higher  centerline  values  (Uc,  Tc)  at  X/D  =  40. 


The  planar  jet  profiles  were  evaluated  at  X/H  =  35  and  are  shown  in  Figures  25  and  26. 
Again,  the  CFD  simulations  display  larger  centerline  values  (Uc,  Tc)  due  to  excessive  virtual 
origins  (see  Figures  17  and  18).  Even  though  the  planar  jet  simulation  predicted  spreading  and 
centerline  decay  rates  (K’s)  closer  to  experimental  data  than  the  round  jet  simulation,  the 
differences  in  the  actual  profiles  at  X/H  =  35  was  generally  worse  for  the  planar  jet  due  to  its 
excessive  virtual  origins.  At  different  distances  from  both  jet  exits,  the  actual  results  will  vary 
due  to  the  offsetting  errors  in  the  virtual  origins  and  spreading/decay  rates. 
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Figure  23.  Actual  Round  Jet  Velocity  Profiles  at  X=1.016m  (X/D=40) 


Figure  24.  Actual  Round  Jet  Temperature  Profiles  at  X=1.016m  (X/D=40) 
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Figure  25.  Actual  Planar  Jet  Velocity  Profiles  at  X=0.455m  (X/H=35) 


Figure  26.  Actual  Planar  Jet  Temperature  Profiles  at  X=0.455m  (X/H=35) 
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VI.  SENSITIVITY  ANALYSIS 

To  evaluate  the  sensitivity  of  the  CFD  models  and  possibly  produce  more  accurate 
results,  various  parameters  were  changed  to  observe  their  resultant  effects.  This  baseline 
knowledge  will  be  used  in  further  research  to  help  predict  jet-type  flows  more  accurately. 

A.  JET  TURBULENCE  INTENSITY 

The  jet  turbulence  intensity  was  varied  from  0.002  to  0.01  on  the  round  jet  (nominal  = 
0.005)  and  from  0.001  to  0.01  on  the  planar  jet  (nominal  =  0.002).  Table  9  and  Figures  27  to  34 
in  Appendix  D  (Section  A)  display  the  results.  The  round  jets  virtual  origins  (C’s)  decrease 
slightly  with  increasing  jet  turbulence  intensity  while  the  multiplication  constants  (K’s)  remain 
relatively  constant.  The  planar  jet  displays  similar  characteristics,  but  deviated  in  different 
directions.  The  smaller  differences  may  be  due  to  slight  variations  in  the  regression  evaluation 
and  not  a  result  of  jet  turbulence  intensity  changes.  Overall,  the  effects  were  minimal  (<5%) 
within  the  given  turbulence  range. 

B.  TURBULENT  PRANDTL  NUMBER 

The  turbulent  Prandtl  Number  (PrO  had  one  of  the  largest  effects  (-15%)  on  the  jet 
predictions.  The  effects  of  changing  the  Prt  by  ±0.2  are  shown  in  Table  10  and  Figures  35  to  42 
in  Appendix  D  (Section  B).  For  both  jets,  the  Prt  had  no  effect  on  the  velocity  characteristics  as 
expected,  but  markedly  changed  the  temperature  features.  When  the  Prt  was  lowered,  the  CFD 
code  correctly  produced  wider  temperature  profiles  (increased  spreading)  and  larger  centerline 
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temperatures  decay  rates  as  more  emphasis  was  artificially  placed  on  the  heat  transport 
mechanism  relative  to  the  momentum  transport  mechanism. 

In  general,  the  Prt  can  be  adjusted  to  produce  the  desired  temperature  characteristics  in 
relation  to  the  velocity  profiles.  For  both  Jet  simulations,  the  recommended  free  jet  Prt’s  (0.5- 
Planar,  0.7-Round)  produced  temperature  profiles  that  were  approximately  same  magnitude 
greater  than  the  velocities  profiles  as  reported  by  experimental  data.  However,  for  the  round  jet 
especially,  the  predicted  velocity  profiles  were  excessively  wide,  therefore  the  resultant 
temperature  profiles  were  also  too  broad. 

C.  JET  LENGTH  SCALE 

Varying  the  characteristic  length  (L)  artificially  changed  the  jet  turbulent  length  scale  (£) 
as  described  in  Appendix  B.  The  characteristic  length  was  varied  from  0.001m  to  0.1m  for  both 
jets  to  investigate  the  resultant  effect.  Table  1 1  and  Figures  43  to  50  in  Appendix  D  (Section  C) 
display  the  results.  As  shown,  changes  in  the  results  were  negligible  (<1%)  and  the  original 
length  scales  used  in  each  simulation  were  determined  to  be  adequate. 

D.  JET  VELOCITY 

The  jet  velocity  was  varied  from  20  to  60  m/s  on  the  round  jet  (nominal  =  40m/s)  and 
between  15  and  55  m/s  on  the  planar  jet  (nominal  =  35m/s).  These  ranges  were  limited  to  typical 
experimental  values  in  order  to  identify  any  trends  at  these  lower  velocities.  The  results  are 
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shown  in  Table  12  and  Figures  51  through  58  in  Appendix  D  (Section  D).  The  change  in 
velocity  did  not  appreciably  effect  (<1%)  either  jet  prediction.  Since  no  boundary  layers 
alterations  were  simulated  at  the  jet  nozzle  (top  hat  velocity  profile),  the  velocity  variations  did 
not  change  the  solution  results. 

E.  JET  TEMPERATURE 

The  jet  exit  temperatures  for  both  models  were  elevated  to  500K  degrees  (200K  above 
ambient  surroundings)  to  show  the  effects  of  a  larger  temperature  difference.  By  limiting 
differential  temperatures  to  less  than  200K,  errors  associated  with  using  the  “constant  specific 
heat  capacity”  option  within  the  simulation  were  minimized.  As  mentioned  before,  the  “constant 
specific  heat  capacity”  option  greatly  reduced  solution  convergence  problems  and  is  consistent 
with  low  working  differential  temperatures. 

The  results  are  displayed  in  Table  13  and  Figures  59  to  66  in  Appendix  D  (Section  E). 
As  indicated  the  temperature  difference  effects  on  the  velocity  and  temperature  profiles  (and 
spreading  rate)  are  negligible  (<1%)  within  the  simulated  temperature  range.  However,  the 
velocity  and  temperature  centerline  decay  rates  (K2a,  K2O  increased  markedly  (~5%).  The 
centerline  decay  characteristics  are  consistent  with  experimental  data  as  described  in  Chapter  HI, 
Section  C.2,  however  the  simulated  spreading  rates  (Kiu,  K;f)  do  not  exhibit  the  coherent 
stracture  effect,  as  expected. 
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F. 


GRID  REFINEMENT 


The  individual  cell  sizes  for  the  round  jet  model  was  reduced  by  two-thirds  (three-fourths 
for  planar  jet)  to  investigate  grid  dependence  of  the  solution.  Since  the  computational  ability  was 
limited,  overall  grid  length  was  reduced  (by  half)  to  keep  the  total  number  of  cells  less  than 
300,000.  All  other  aspects  of  the  models  were  kept  the  same.  Figures  67  to  74  in  Appendix  D 
(Section  F)  illustrate  the  results.  The  plots  are  shown  differently  due  to  the  inability  of  the 
shorter  refined  model  to  reach  similarity.  Therefore,  the  similarity  equations  could  not  be  used 
to  evaluate  the  jet  characteristics  and  the  predictions  had  to  be  displayed  as  shown. 

The  differences  between  the  nominal  and  refined  models  were  found  to  be  small  (<5%), 
therefore  grid  resolution  was  deemed  satisfactory  for  both  jets.  Several  additional  grid 
configurations  were  attempted  without  successful  solution  convergence.  However,  future  grid 
refinement  in  the  jet  shear  layer  and  other  high  gradient  areas  may  be  able  to  produce  more 
significant  changes  in  the  solution  results.  Time  constraints  limited  further  research  in  this  area. 
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VII.  SUMMARY  AND  RECOMMENDATIONS 


A  critical  review  of  free-jet  experimental  data  was  summarized  and  compared  to  CFD 
simulations.  The  CFDRC  flow  code,  using  the  k-e  turbulence  model,  predicted  round 
axis5rmmetric  jets  poorly  while  marginally  estimating  the  two-dimensional  planar  jet  flow.  By 
not  being  able  to  vary  some  of  the  turbulence  model  constants,  the  code  over  predicted  the 
velocity  spreading  and  centerline  decay  of  the  round  jet  by  approximately  25  percent.  Accuracy 
of  the  planar  jet  simulations  was  generally  better  within  an  approximately  15  overall  percent 
deviation,  however  larger  centerline  temperature  decay  errors  existed. 

Sensitivity  analysis  on  both  jet  simulations  indicated  that  changes  in  jet  exit  velocity 
(Ue),  turbulence  intensity  (Tu),  and  turbulent  length  scale  (L)  had  little  or  no  effect  on  the 
solution  results.  The  jet  exit  temperamre  (Te)  and  turbulent  Prandtl  Number  (PrO,  conversely, 
produced  significant  changes  in  the  results.  The  recommended  free-jet  Prandtl  Numbers 
(Chapter  HI,  Section  E.2)  resulted  in  correct  temperature  profiles  (spreading  rates)  when 
compared  to  the  velocity  profiles.  The  change  in  jet  exit  temperature  also  created  correct 
centerline  decay  variations,  but  failed  to  produce  appropriate  spreading  rate  changes  as  reported 
in  experiments. 

Even  though  all  attempts  to  increase  the  accuracy  of  the  simulations  failed,  the  sensitivity 
analysis  and  insights  reported  in  this  study  should  be  extremely  useful  to  follow-on  researchers 
utilizing  CFDRC’ s  CFD-ACE-i-  code.  Recommendations  for  further  work  include  evaluating  the 
alternative  RNG  k-e  turbulence  model  contained  within  the  program.  This  turbulence  model  is 
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an  updated  variation  of  the  k-z  model  that  uses  a  renormalization  group  approach  to 
systematically  remove  the  smallest  scales  of  motion.  The  RNG  k-z  turbulence  model  was  tried 
several  times  without  obtaining  proper  solution  convergence,  most  probably  due  to  grid 
size/refmement  limitations  (<300,000  cells).  Time  constraints  prevented  further  in-depth 
troubleshooting  attempts 

More  realistic  jet  exit  profiles  should  also  be  simulated  and  compared  to  experimental 
data.  By  using  Fortran  input  subroutines  within  the  code,  the  user  can  import  realistic  jet 
velocity  and  temperature  exit  profiles.  Since  jet  nozzle  boundary  layers  have  a  significant  effect 
on  real  jet  behavior,  their  effects  and  sensitivity  on  the  simulated  results  should  also  be 
investigated.  In  addition,  turbulence  quantities  (^,8)  of  the  initial  conditions  should  be  varied  to 
investigate  their  possible  effects. 
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APPENDIX  A.  THE  CFDRC  PROGRAM 


The  CFD  Research  Corporation  (CFDRC)  provides  a  variety  of  tools  for  the  simulation 
and  analysis  of  fluid  flow  and  associated  physics  for  an  assortment  of  industrial  applications. 
CFDRC  has  specifically  developed  their  software  so  that  the  average  engineering  professional 
can  easily  manipulate  several  multi-disciplinary  engineering  project  simulations. 

Typical  numerical  simulations  involve  three  distinct  process  steps: 

1)  The  volume  of  interest  (or  solution  space)  must  be  divided  into  discrete  control  volumes  or 
cells. 

2)  One  must  then  define  the  boundary  conditions,  initial  conditions,  and  the  required  equations 
to  be  solved  at  each  cell.  In  addition,  the  numerical  technique  used  to  solve  the  required 
equations  must  also  be  defined. 

3)  Finally,  after  the  solution  has  been  calculated,  the  information  needed  must  be  extracted  from 
the  large  volume  of  data  generated  in  the  solution  process. 

CFDRC  provides  software  modules  and  packages  to  address  each  of  these  steps  for  a 
total  overall  solution.  This  program,  entitled  CFD-ACE+  (Version  6.2),  uses  three  separate,  yet 
interactive  codes  to  impart  an  “all  in  one”  seamless  commercial  code  flow  solver.  CFD-GEOM 
provides  interactive  geometry  modeling  and  grid  generation  capabilities.  Similar  to  computer 
aided  design  (CAD),  geometry  modeling  is  the  process  of  creating  a  computer  model  of  the 
geometry  that  makes  up  a  problem  domain.  CFD-GEOM  offers  a  NURBS-based  geometry 
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engine  with  a  variety  of  geometric  construction  tools.  It  also  has  the  ability  to  import  and  export 
IGES  data  from  most  major  industry  CAD  programs.  Grid  generation  is  the  process  of 
discretizing  the  problem  domain  with  individual  cells  over  which  the  flow  equations  are 
integrated.  CFD-GEOM  provides  for  the  production  of  two  classes  of  cells:  structured  and 
unstructured. 

CFD-ACE  (U)  is  an  unstructured,  polyhedral  cell  flow  solver.  It  is  also  integrated  with  a 
wide  variety  of  physics  modules  making  it  the  core  of  a  multi-disciplinary  analysis  environment. 
Inputs  are  specified  for  CFD-ACE  (U)  using  CFD-GUI,  an  advanced  graphical  user  interface  that 
allows  complete  specification  of  the  multi-physics  problem. 

CFD-ACE  (U)  employs  a  cell-centered  control  volume  solution  approach.  This  approach 
implies  that  the  discrete  equations  are  formulated  by  evaluating  and  integrating  the  fluxes  across 
the  faces  that  surround  each  control  volume.  In  addition,  CFD-ACE  (U)  uses  a  pressure-based 
methodology  in  which  pressure  becomes  one  of  the  dependent  variables  evaluated  at  each  cell. 

The  CFD-ACE  (U)  unstructured  flow  solver  can  simulate  a  wide  variety  of  flow  regimes 
and  phenomenon  as  listed  below. 

•  Internal  or  External  Flow 

•  Laminar  or  Turbulent  Flow:  (involves  the  solution  of  one  or  more  additional  equations) 

•  Incompressible  or  Compressible  Flow 
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Heat  Transfer 


•  Mixing  and  Reaction;  Flows  involving  multiple  gases  that  mix 

and  react  can  be  modeled. 

•  Steady-  state  or  Transient 

•  Several  others  [see  CFDRC  User  Manuals  (2000)] 

CFD-ACE  (U)  also  offers  a  variety  of  differencing  schemes  (1® -order  upwind  to  high¬ 
speed  S'^-order)  to  compute  the  variables  within  the  domain.  The  type  of  differencing  scheme 
determines  how  the  cell  face  values  are  calculated  from  adjacent  cells.  These  schemes  can  be 
independently  selected  for  each  variable  to  be  solved;  however,  for  this  investigation  a  U'-order 
upwind  scheme  was  chosen  for  all  variables  in  the  jet  models  to  simplify  convergence 
operations. 

CFD-ACE  (U)  employs  an  iterative  solution  technique  in  which  the  assembled  equations 
for  each  dependent  variable  are  solved  sequentially  and  repeatedly,  with  the  goal  of  improving 
the  overall  solution  with  each  iteration.  This  solution  convergence  reduces  the  variable  value 
changes  within  each  control  volume  until  acceptably  small  values  are  obtained  and  the  final 
overall  solution  is  determined.  The  nonlinear  coupled  nature  of  the  Navier-Stokes  (and  other 
relevant  equations)  makes  it  necessary  to  restrain  or  under-relax  the  iteration-to-iteration  changes 
of  each  variable  in  order  to  obtain  a  stable  convergence  of  the  solution  procedure.  Under¬ 
relaxation  constrains  the  amount  that  each  variable  can  change  from  one  iteration  to  the  next. 
The  dependent  variables  (u,  v,  w,  and  k)  are  modified  by  using  an  Inertial  Factor.  However,  a 
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linear  under-relaxation  technique  is  applied  to  the  auxiliary  variables  (p,  T,  }i,  p).  The  code 

allows  the  user  to  change  these  under-relaxation  values  within  the  program  to  ensure  stable 
solution  convergence. 

CFD-ACE  (U)  generates  a  wide  variety  of  graphical  and  printed  outputs.  This  output 
includes  graphics  files  with  information  at  each  cell  for  use  in  CFD-VIEW  and  the  printed 
output,  which  allows  the  user  to  monitor  several  variables  at  a  fixed  location  in  the  problem 
domain.  Integrated  quantities  such  as  mass  flow-rates,  heat  transfer  rates,  and  pressure  forces 
can  be  written  to  the  text  output  file.  In  addition,  the  change  in  the  solution  from  iteration  to 
iteration,  or  the  residual,  can  be  graphically  monitored  to  assess  convergence. 

One  of  the  larger  challenges  in  computational  modeling  is  the  management  of  the  large 
volume  of  data  generated  with  each  simulation.  These  data  must  be  reduced  to  extract  useful 
information,  which  can  be  applied  to  practical  problems.  To  aid  in  the  data  reduction  process, 
the  CFD-ACE+  suite  includes  a  3-D  graphical  post-processor  called  CFD-VIEW. 

CFD-VIEW  contains  a  variety  of  tools  to  visualize  and  extract  data  from  complex  3-D 
data  sets.  Various  types  of  visual  surfaces  can  be  generated  such  as  constant  computational 
plane  surfaces,  cutting  plane  surfaces,  iso-value  surfaces,  and  unstructured  surfaces.  On  each 
surface,  contour  levels  of  any  of  the  variables  in  the  data  set  can  be  displayed  as  value-colored 
lines,  flooded  contour  levels,  and  continuously  shaded  value-colored  surfaces.  Vector  fields  can 
also  be  displayed  on  these  planes  using  arrows  that  indicate  magnitude  and  direction. 
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Other  miscellaneous  features  of  CFD-VIEW  include  a  point  probe  that  allows  one  to 
extract  data  from  any  point  in  the  flow  field,  a  line  probe  used  to  make  X-Y  plots,  and  a 
streamline  tracer  to  create  particle  traces  using  any  vector  field  in  the  model.  CFD-VIEW  also 
contains  an  expression  calculator  that  allows  one  to  perform  many  of  the  tasks  necessary  to 
reduce  a  computational  data  set  into  usable  information.  The  calculator  has  a  variety  of  built-in 
functions  that  enable  cell-wise  mathematical  operations  and  derivative  mathematical  options  that 
depend  on  more  than  one  cell.  Results  of  calculations  can  be  visualized  in  the  same  way  as  any 
other  data. 

Overall,  the  three  interactive  codes  of  CFD-ACE+  {CFD-GEOM,  CFD-GUI  /  CFD-ACE 
(U),  and  CFD-VIEW}  work  in  harmony  to  step  the  user  through  the  numerical  simulation 
process.  The  reader  is  referred  to  the  CFDRC  User  Manuals  (2000)  for  more  comprehensive 
information. 
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APPENDIX  B.  CFD-ACE  (U)  STANDARD  k-e  TURBULENCE  MODEL 

The  non-linearity  of  the  Navier-Stokes  equations,  coupled  with  the  complexity  of  the 
boundary  conditions,  makes  it  impossible  to  obtain  analytical  solutions  for  all  but  a  limited 
number  of  flows  of  engineering  interest.  Hence  one  is  forced  to  resort  to  approximate  or 
numerical  methods.  Even  though  a  wide  variety  of  numerical  techniques  can  be  applied  to  solve 
the  Navier-Stokes  equations  for  laminar  flows.  Direct  Numerical  Simulation  (DNS)  of  turbulent 
flows  is  feasible  only  at  very  low  Reynolds  numbers.  Turbulent  flows  are  inherently  unsteady 
and  they  contain  a  wide  range  of  time  and  length  scales,  and  resolution  of  these  scales  requires 
very  short  time  steps  and  fine  grids.  The  CPU  and  memory  requirements  are  too  large  for  most 
present  day  computers. 

As  most  engineering  applications  only  require  time-mean  quantities,  the  Navier-  Stokes 
equations  are  usually  averaged  over  time  or  ensemble  of  statistically  equivalent  flows  to  yield 
averaged  equations,  hi  the  averaging  process,  a  flow  quantity  is  decomposed  into  mean  and 
fluctuating  parts.  Reynolds  (or  time)  averaging  and  Favre  (or  density)  averaging  are  two  of  the 
more  popular  techniques  generally  used. 

For  turbulent  flows,  CFD-ACE  (U)  exclusively  utilizes  the  Favre-Averaged  Navier- 
Stokes  (FANS)  equations  to  solve  for  momentum  and  energy.  Although  the  FANS  equations 
contain  less  information  than  the  full  NS  equations,  they  do  contain  additional  unknown 
Reynolds  stresses.  These  correlations  between  the  fluctuating  components  arise  in  the  averaging 
process,  and  are  additionally  modeled  to  achieve  closure  of  the  FANS  equations.  All  the 
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turbulence  models  available  in  CFD-ACE  (U)  employ  the  generalized  Boussinesq  eddy  viscosity 
concept  in  which  the  Reynolds  stresses  are  treated  as  a  linear  function  of  the  mean  strain  rate. 
The  Generalized  Transport  Equation  given  below  indicates  the  common  form  for  all  of  the 
FANS  equations. 


—{p<l))  +  —  (pu.<p)= 
Ot  dX; 


— r  -^+5 


(38) 


{Transient}  (Convection)  (Diffusion)  (Source) 


The  symbol  tprazy  represent  any  of  the  velocity  components,  enthalpy,  or  other  scalar 
flow  variables.  In  the  preceding  equation,  Fg^is  the  effective  diffusivity  and  is  modeled  as: 

r  -JL+JL.  (39) 

^^“Pr  Pr, 

is  a  generalized  source  term,  which  represents  the  mechanisms  for  the  generation  and 
destruction  of  <}>.  In  addition,  any  terms  that  cannot  be  conveniently  expressed  as  convection  or 
diffusion  (e.  g.,  the  pressure  gradient  term  in  the  momentum  equations)  can,  in  general,  be 
lumped  into  the  source  term  (S^). 


The  standard  k-z  model,  employed  by  CFD-ACE  (U),  is  based  on  Launder  /  Spalding 
(1974).  The  two-equation  model,  in  a  FANS  generalized  form,  governs  the  transport  of  turbulent 
kinetic  energy  (k)  and  its  dissipation  rate  (e).  The  square  root  of  k  is  taken  to  be  the  velocity 
scale,  while  the  length  scale  {t)  is  modeled  as: 


1  = 


e 


(40) 
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The  expression  for  eddy  kinematic  viscosity  is: 


'  p  « 


Therefore,  the  modeled  equations  for  k  and  e  are: 


i  LV 
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with  the  production  (P)  defined  as: 

f du;  ^Uj  2  du 


P  =  v, 
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(44) 


The  five  constants  used  in  the  fc-f  model  are  listed  below.  All  of  them  are  internal 
to  the  code  and  can  not  be  changed  for  axisymmetrical  flows  as  suggested  by  Launder/Spalding 


C^=0.09 
C„=1.44 
C,2=1.92 
a,  =1.0 
<7^  =1.3 

(1974).  Therefore,  without  being  able  to  vary  and  Qi  across  the  jet  profile,  the  overall 
prediction  accuracy  for  axisymmetrical  jets  will  be  poor. 
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The  standard  k-z  model  is  a  high  Reynolds  Number  model  and  is  not  intended  to  be  used 
in  the  near-wall  regions  where  viscous  effects  dominate  the  effects  of  turbulence.  Instead,  “wall 
functions”  are  used  in  cells  adjacent  to  walls.  The  k  and  8  transport  equations  are  not 
numerically  integrated  in  these  cells.  Instead,  semi-empirical  expressions  are  used  to  relate  k, 
s,  and  the  friction  velocity  (w).  These  expressions  are  obtained  from  analysis  of  the  momentum 
and  turbulence  equations  for  a  flat  plate  boundary  layer,  assuming  a  logarithmic  velocity  profile. 


(45) 

cfky^ 

(46) 

The  friction  velocity  (u)  is  defined  as: 

i~p  (47) 


where  xw  is  the  shear  stress  at  the  wall  and  is  obtained  by  assuming  the  velocity  profile  between 
the  wall  and  the  next-to-wall  grid  points  obeys  the  following  “law  of  the  wall”: 


For  <  11.5 

(48) 

For  />11.5 

+ 

11 

•f 

(49) 
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with  the  dimensionless  distances;'^  and  defined  as; 

(50) 

M 

(51) 

w 

The  wall  shear  stress  (xw)  is  calculated  iteratively  from  the  known  values  of:F  and  «  in 
the  first  cell.  The  constants  appearing  in  Equations  46  and  49  are  experimentally  determined  to 
be  E  =9.  0  and  K  =  0.4  (von-Karman  constant).  Because  the  semi-empiricdl  relations  for  k  and 
8  in  the  first  cell  assume  a  logarithmic  velocity  profile,  the  turbulence  wall  functions  are  strictly 
valid  only  if  the  center  of  the  cell  nearest  the  wall  is  inside  the  logarithmic  boundary  layer  (30  < 
/  <  150  suggested  range).  The  wall  shear  stress,  evaluated  using  Equations  48  and  49,  is  used  to 
calculate  the  boundary  condition  for  the  velocity  components  parallel  to  the  wall. 

For  specifying  tuibulence  characteristics  at  inlet  and  outlet  boundaries,  three  quantities 
can  be  utilized;  turbulent  kinetic  energy  (A:),  dissipation  rate  (e),  and  length  scale  (f).  The  code 
requires  that  only  two  of  the  three  quantities  be  specified.  It  is  usually  more  convenient,  to 
provide  a  length  scale  instead  of  a  dissipation  rate  value.  For  example,  the  length  scale  is  usually 
the  inlet  diameter  or  height.  In  this  case,  the  program  will  internally  calculate  the  boundary 
value  of  8  from  the  given  values  of  k  and  i  as  shoAvn  in  Equation  52,  where  the  constant,  still 
has  the  value  0.09. 
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This  formula  assumes  that  the  characteristic  length  scale  for  the  turbulent  eddies  is 
approximately  0.3%  of  the  size  of  the  inlet  dimension  [L  =  D  (or  H)  in  our  case].  The  eddy- 
viscosity  (ji,)  at  the  boundary  is  then  internally  calculated  using  the  given  k  and  the  calculated  (or 
given)  e  as: 


= 


C,pk^ 

e 


(53) 


After  applying  the  Favre  averaging  procedure,  the  total  enthalpy  energy  equation 
becomes: 


K 


ar 

a^. 


a 

_i _ 

at 

'dx., 

(54) 


where  the  total  (stagnation)  enthalpy  is  defined  as: 


Similar  to  the  momentum  equation  analysis,  the  viscous  stress  tensor  (•?;)  for  Newtonian 
fluids  can  be  related  to  the  velocity  gradients  by  the  following  equation. 


=/^ 


2  (du  ^ 

-+  •' 


duj 
dX;  dX; 


‘ ) 


— 

3  ' 


V 


dx^ 


f'lj 


where  4  is  the  Kroneker  delta  function: 


for  i  ^  j 
for  i  =  j 


(56) 
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APPENDIX  C.  SAMPLE  SPREADSHEET 
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0.2320681  300.1876 


1.096818  6.99516  303.9301  43.18181102  8.872002082  7.833393562 

1.106364  5.913419  303.8712  43.83837796  6.784228951  7.749636028 

1.119909  5.84544  303.8139  44.09090661  6.842892922  7.886983974 

1.131455  5.772453  303.7627  44.64647244  6.929414683  7.972998113 

1.143  5.707756  303.7317  46  7.007968996  8.039231449 


APPENDIX  D.  SENSITIVITY  ANALYSIS  nCURES  AND  TABLES 


A.  JET  TURBULENCE  INTENSITY 


Table  9.  CFD  Simulation  Predictions  with  Varying  Jet  Turbulence  Intensity 


Jet 

Turb. 

Intensity 

Vel. 

Spread 

Rate 

Geom. 

Virt 

Orig. 

C/L 

Vel. 

Decay 

Kin. 

Virt 

Orig. 

Temp. 

Spread 

Rate 

Spread 

Virt. 

Orig. 

CTL 

Temp. 

Decay 

Decay 

Virt 

Orig. 

{Tu} 

{K,J 

(C.J 

{K2„} 

{C2„} 

{KnI 

{C,.} 

IKz.) 

fC2tl 

1  AXISYMMETRICAL  ROUND  JET  | 

Target 

Values 

0.085- 

0.095 

-1.0  to 
4.0 

0.160- 

0.172 

-1.0  to 
4.0 

0.098- 

0.110 

0.0  to 
2.0 

0.185- 

0.194 

0.0  to 
2.0 

0.002 

0.111 

9.4 

0.200 

10.2 

0.246 

12.5 

Him 

9.3  n 

0.201 

10.0 

0.134 

11.1 

0.247 

12.3 

0.115 

8.8  ^ 

0.203 

0.249 

11.4 

1  TWO-DIMENSIONAL  PLANAR  JET  | 

-6.0  to 
2.0 

0.155- 

0.220 

-2.0  to 
4.0 

0.125- 

0.140 

-5.0  to 
0.0 

0.001 

3.6 

0.152 

5.1 

0.145 

8.4 

0.207 

9.1 

V7iTiTi7>!1llll 

3.7 

0.154 

5.2  1 

0.149 

8.6 

0.209 

9.1 

0.01 

0.092 

3.3 

0.154 

4.8 

0.149 

8.4 

0.208 

8.5 

87 


(Te-Tinf)/(dTc-Tinf)  ^  Ue/Uc 


X/D 


Figure  28.  Round  Jet  Centerline  Temp.  Decay  with  Varying  Jet  Turb.  Intensity 
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Eta=FV(X-C1uD) 


Figure  29.  Round  Jet  Velocity  Profile  with  Varying  Jet  Turbulence  Intensity 


Eta=IVpC-C1tD) 

Figure  30.  Round  Jet  Temp.  Profile  with  Varying  Jet  Turbulence  Intensity 
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(Te-Tlnf)'^2/(Tc-Tlnf)''2  ^  (Ue/Uc)'^2 


X/H 


Figure  33.  Planar  Jet  Velocity  Profile  with  Varying  Jet  Turbulence  Intensity 


Figure  34.  Planar  Jet  Temp.  Profile  with  Varying  Jet  Turbulence  Intensity 
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B 


TURBULENT  PRANDTL  NUMBER 


Table  10.  CFD  Simulation  Predictions  with  Varying  Prandtl  Number 


Jet 

Vel. 

Geom. 

C/L 

Kin. 

Temp. 

Spread 

C/L 

Decay 

Turb. 

Spread 

Virt. 

Vel. 

Virt. 

Spread 

Virt 

Temp. 

Virt. 

Pr. 

Rate 

Orig. 

Decay 

Orig. 

Rate 

Orig. 

Decay 

Orig. 

{Prt} 

{K,J 

{C,„} 

{KaJ 

{Ca„} 

IKu) 

{Cul 

{Ka.} 

{Ca.) 

AXISYMMETRICAL  R< 

OUND  JET 

Target 

0.096- 

-1.0  to 

0.155- 

-1.0  to 

0.125- 

0.0  to 

0.250- 

0.0  to 

Values 

0.120 

4.0 

0.220 

4.0 

0.140 

2.0 

0.280 

2.0 

0.5 

0.112 

9.3 

0.201 

10.0 

0.159 

11.4 

0.286 

12.7 

(0.7)„„„, 

0.112 

9.3 

0.201 

10.0 

0.134 

11.1 

0.247 

12.3 

0.9 

0.111 

9.2 

0.199 

10.0 

0.112 

9.4 

0.214 

11.5 

TWO-DIMENSIONAL  PLANAR  JET 

Target 

0.096- 

-6.0  to 

0.155- 

-2.0  to 

0.125- 

-5.0  to 

0.250- 

-5.0  to 

Values 

0.120 

2.0 

0.220 

4.0 

0.140 

0.0 

0.280 

0.0 

0.3 

0.091 

3.6 

0.152 

5.1 

0.173 

5.6 

0.240 

8.4 

(0.5)nom 

0.092 

3.7 

0.154 

5.2 

0.149 

8.6 

0.209 

9.1 

0.7 

0.091 

3.6 

0.152 

5.1 

0.118 

7.2 

0.177 

8.1 
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Figure  35.  Round  Jet  Centerline  Vel.  Decay  with  Varying  Prandtl  Number 


Figure  36.  Round  Jet  Centerline  Temp.  Decay  with  Varying  Prandtl  Number 
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Eta=R/(X-C1uD) 


Figure  37.  Round  Jet  Velocity  Profile  with  Varying  Prandtl  Number 


Eta=R/(X-C1tD) 

Figure  38.  Round  Jet  Temp.  Profile  with  Varying  Prandtl  Number 
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Figure  39.  Planar  Jet  Centerline  Vel.  Decay  with  Vaiying  Prandtl  Number 


Figure  40. 


Planar  Jet  Centerline  Temp.  Decay  with  Varying  Prandtl  Number 


T-Tinf)/(Tc-Tinf 


Ba=Y/(X-C1tH) 

Figure  42.  Planar  Jet  Temp.  Profile  with  Varying  Prandtl  Number 
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C.  JET  LENGTH  SCALE 


Table  11.  CFD  Simulation  Predictions  with  Varying  Jet  Length  Scale 


Jet 

Length 

Scale 

Vel. 

Spread 

Rate 

Geom. 

Virt. 

Orig. 

C/L 

Vel. 

Decay 

WIJH 

mm 

Temp. 

Spread 

Rate 

Spread 

Virt. 

Orig. 

C/L 

Temp. 

Decay 

Decay 

Virt. 

Orig. 

{Clu}_ 

{K2u} 

IHi 

(K..> 

...{£»L, 

{K2.} 

AXISYMMETl 

HCAL  R( 

OUNDJET 

Target 

Values 

0.085- 

0.095 

-1.0  to 
4.0 

0.160- 

0.172 

-1.0  to 
4.0 

0.098- 

0.110 

0.0  to 
2.0 

0.185- 

0.194 

0.0  to 
2.0 

0.1 

0.111 

9.1 

0.199 

10.0 

0.133 

11.0 

0.245 

12.3 

0.112 

9.3 

0.201 

10.0 

0.134 

11.1 

0.247 

12.3 

0.001 

o.iio" 

9.2 

0.199 

0.132 

0.244 

12.5 

1  TWO-DIMENSIONAL  PLANAR  JET 

Target 

Values 

0.096- 

0.120 

-6.0  to 
2.0 

0.155- 

0.220 

-2.0  to 
4.0 

0.125- 

0.140 

-5.0  to 
0.0 

0.250- 

0.280 

-5.0  to 
0.0 

0.1 

0.091 

3.6 

0.152 

5.1 

0.146 

8.4 

0.207 

9.0 

0.092 

3.7 

5.2 

0.149 

8.6 

0.209 

9.1 

mm\m 

0.091 

3.6 

■rtteVJI 

5.1 

0.145 

8.4 

0.207 

9.1 
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Figure  43.  Round  Jet  Centerline  Vel.  Decay  with  Varying  Jet  Length  Scale 


Figure  44.  Round  Jet  Centerline  Temp.  Decay  with  Varying  Jet  Length  Scale 
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Figure  45.  Round  Jet  Velocity  Profile  with  Varying  Jet  Length  Scale 


Figure  46. 


Round  Jet  Temp.  Profile  with  Varying  Jet  Length  Scale 


(Te-Tlnf)^2/(Tc-Tlnf)^2  'S  (Ue/Uc)'^2 


X/H 

47.  Planar  Jet  CenterBne  Vel.  Decay  with  Varying  Jet  Length  Scale 


X/H 

Figure  48.  Planar  Jet  Centerline  Temp.  Decay  with  Varying  Jet  Length  Scale 


»!« 


Figure  49.  Planar  Jet  Velocity  Profile  with  Varying  Jet  Length  Scale 


Figure  50.  Planar  Jet  Temp.  Profile  with  Varying  Jet  Length  Scale 
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D.  JET  VELOCITY 


Table  12.  CFD  Simulation  Predictions  with  Varying  Jet  Velocity 


Jet 

Vel. 

Geom. 

C/L 

Kin. 

Temp. 

Spread 

C/L 

Decay 

Exit 

Spread 

Virt. 

Vel. 

Virt. 

Spread 

Virt 

Temp. 

Virt. 

Vel 

Rate 

Orig. 

Decay 

Orig. 

Rate 

Orig. 

Decay 

Orig. 

{Ue} 

{K,„} 

(C,„} 

{K^J 

{C2„} 

{K,.} 

{C,J 

{K«} 

{C2.} 

AXISYMMETRICAL  R( 

OUNDJET 

Target 

0.085- 

-1.0  to 

0.160- 

-1.0  to 

0.098- 

0.0  to 

0.185- 

0.0  to 

Values 

0.095 

4.0 

0.172 

4.0 

0.110 

2.0 

0.194 

2.0 

20 

0.111 

9.1 

0.194 

9.5 

0.130 

10.8 

0.250 

12.4 

(40)„om 

0.112 

9.3 

0.201 

10.0 

0.134 

11.1 

0.247 

12.3 

60 

0.113 

9.4 

0.203 

10.2 

0.136 

11.1 

0.240 

12.2 

1  TWO-DIMENSIONAL  PLANAR  JET 

Target 

0.096- 

-6.0  to 

0.155- 

-2.0  to 

0.125- 

-5.0  to 

0.250- 

-5.0  to 

Values 

0.120 

2.0 

0.220 

4.0 

0.140 

0.0 

0.280 

0.0 

15 

0.091 

3.5 

0.153 

5.1 

0.146 

8.3 

0.212 

9.1 

(35)„„m 

0.092 

3.7 

0.154 

5.2 

0.149 

8.6 

0.209 

9.1 

55 

0.091 

3.6 

0.152 

5.1 

0.146 

8.4 

0.198 

8.8 
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X/D 


Figure  51.  Round  Jet  Centerline  Vel.  Decay  with  Varying  Jet  Velocity 
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(T-Tinf)/(Tc-Tinf)  U/Uc 


Eta=R/(X-C1uD) 


Figure  53.  Round  Jet  Velocity  Profile  with  Varying  Jet  Velocity 


Figure  54.  Round  Jet  Temp.  Profile  with  Varying  Jet  Velocity 


104 


(Te-Tinf)''2/(Tc-Tlnl)''2  S  (Ue/Uc)''2 
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Figure  57.  Planar  Jet  Velocity  Profile  with  Vaiying  Jet  Velocity 


Figure  58.  Planar  Jet  Temp.  Profile  with  Varying  Jet  Velocity 
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JET  TEMPERATURE 


Table  13.  CFD  Simulation  Predictions  with  Jet  Temperature 


Jet 

Vel. 

Geom. 

C/L 

Kin. 

Temp. 

Spread 

C/L 

Decay 

Exit 

Spread 

Virt 

Vel. 

Virt. 

Spread 

Virt. 

Temp. 

Virt. 

Temp 

Rate 

Orig. 

Decay 

Orig. 

Rate 

Orig. 

Decay 

Orig. 

{Tel 

{Kiul 

{C,„} 

{KiJ 

{CaJ 

{K„} 

{C,.} 

{Kz.} 

{C.J 

AXISYMMETl 

UCAL  ROUND  JET  | 

Target 

0.160- 

-1.0  to 

0.098- 

0.0  to 

0.185- 

0.0  to 

Values 

0.172 

4.0 

0.110 

2.0 

0.194 

2.0 

(330)„„ 

0.112 

9.3 

0.201 

10.0 

0.134 

11.1 

0.247 

12.3 

500 

9.2 

0.244  1 

10.4 

0.133 

10.8 

0.309 

12.6 

1  TWO-DIMENSIONAL  PLANAR  JET 

Target 

0.096- 

-6.0  to 

0.155- 

-2.0  to 

0.125- 

-5.0  to 

0.250- 

-5.0  to 

Values 

0.120 

2.0 

0.220 

4.0 

0.140 

0.0 

0.280 

0.0 

3.7 

0.154 

5.2 

8.6 

0.209 

9.1 

500 

0.093 

4.2 

0.148 

8.4 

10. 

0.301 

10.5 
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(T-Tinf)/(Tc-Tlnf)  U/Uc 


Figure  61.  Round  Jet  Velocity  Profile  with  Varying  Temperature 


Figure  62.  Round  Jet  Temp.  Profile  with  Varying  Temperature 
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X/H 

Figure  63.  Planar  Jet  Centerline  Vel.  Decay  with  Varying  Temperature 


XAH 

Figure  64.  Planar  Jet  Centerline  Temp.  Decay  with  Varying  Temperature 


no 


Figure  65.  Planar  Jet  Velocity  Profile  with  Varying  Temperature 


Figure  66.  Planar  Jet  Temp.  Profile  with  Varying  Temperature 
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F. 


GRID  REFINEMENT 


Figure  67.  Round  Jet  Centeriine  Vel.  Decay  with  Grid  Refinement 


Figure  68.  Round  Jet  Centerline  Temp.  Decay  with  Grid  Refinement 


112 


T  (K)  U  (m/s) 


Figure  69.  Round  Jet  Velocity  Profile  with  Grid  Refinement 


Temperature  Profile  @  X=0.5m(20D) 


Figure  70.  Round  Jet  Temp.  Profile  with  Grid  Refinement 
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■c-Tinf)]''2 


0  5  10  15  20  25  30  35  40 

X/H 

Figure  71.  Planar  Jet  Centerline  Vel.  Decay  with  Grid  Refinement 


0  5  10  15  20  25  30  35  40 

XAH 


Figure  72.  Planar  Jet  Centerline  Temp.  Decay  with  Grid  Refinement 
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Velocity  Profile  @  X=0.2m(15H) 


Figure  73.  Planar  Jet  Velocity  Profile  with  Grid  Refinement 


Figure  74.  Planar  Jet  Temp.  Profile  with  Grid  Refinement 
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